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I .     INTRODUCTION 

Most  engineering  structural  applications  involve  the  use 
of  materials  containing  second  phases  in  the  form  of  discrete 
precipitates  or  distributions  of  particulates.  This  latter 
category  includes  discontinuously  reinforced  metal-matrix 
composites.  The  behavior  of  these  materials  during  processing 
and  the  final  mechanical  properties  depend  upon  the  nature, 
volume  fraction,  size  and  also  the  dispersion  of  the 
particles. 

In  both  conventional  alloys  and  in  metal-matrix  composites 
the  second  phase  particles  may  be  obstacles  to  the  movement  of 
dislocations.  If  a  dislocation  is  prevented  from  completely 
traversing  the  crystal  under  an  applied  stress,  then  the 
crystal  will  become  harder  to  deform.  A  higher  stress  is  then 
required  to  move  the  dislocation  past  the  obstacle  and 
therefore  the  material  must  be  stronger.  Orowan  [Ref.  1] 
found  that  the  shear  stress  required  for  a  dislocation  to 
bypass  a  pair  of  small  particles  may  be  expressed  as: 

where  GM  is  the  shear  modulus  of  the  matrix,  b  is  the  burgers 
vector  and  D  is  the  space  between  particles  in  the  slip  plane. 
From  this  relationship  it  can  be  seen  that  a  smaller  inter- 


particle  spacing  will  result  in  a  larger  stress  to  move  the 
dislocation  past  the  particles. 

Larger  second  phase  particles  will  also  affect  nucleation 
and  grain  growth  during  recrystallization  of  conventional 
alloys.  In  general  it  has  been  found  that  greater  inter- 
particle  spacing  and  larger  particles  will  facilitate  the 
process  of  recrystallization  [Ref.  2].  Close  particle  spacing 
and  small  particles  have  the  opposite  effect.  The  second- 
phase  particles  increase  the  dislocation  density  when  the 
material  is  plastically  deformed,  increasing  the  driving  force 
for  recrystallization.  However,  they  also  impede  the 
rearrangement  of  dislocations  to  form  mobil  high  angle  grain 
boundaries  and  thereby  hinder  the  migration  of  these  grain 
boundaries,  thus  retarding  recrystallization. 

In  research  on  Al  alloys  with  an  AlCu2  particle  size 
greater  than  0.5  /um  it  was  found  that  inter-particle  spacing 
had  the  greatest  effect  on  recrystallization  [Ref.  2].  For 
wide  particle  spacings,  subgrains  formed  by  the  rearrangement 
of  dislocations  and  grew  until  they  reached  a  critical  size, 
at  which  point  they  would  begin  to  migrate.  The  movement  of 
the  subgrains  was  impeded  by  the  particles.  This  allowed  for 
further  nucleation  of  additional  subgrains  and  resulted  in 
accelerated  recrystallization. 

The  deformation  structure  of  a  deformed  alloy  is  also 
dependent  on  the  particle  size.   Large  particles  cause  local 


distortion  of  the  lattice  and  increase  the  rate  of  nucleation. 
Small  particles  increase  the  homogeneity  of  the  dislocation 
distribution  which  slows  down  the  nucleation  rate.  Therefore, 
for  an  arrangement  of  coarse,  widely  dispersed  particles  the 
recrystallization  is  significantly  increased.  On  the  other 
hand,  for  a  small  particle  diameter  and  small  inter-particle 
spacing,  the  dislocation  distribution  is  much  more  dense  and 
homogeneous  which  reduces  the  rate  of  recrystallization. 

Second  phase  particle  distributions  may  be  either  random 
or  non-random.  Random  or  uniform  distribution  infers  that 
local  particle  densities  are  the  same  throughout  the  material. 
Non-random  distributions  are  just  the  opposite  and  may  be  the 
result  of  insufficient  processing  of  the  material.  High  and 
low  density  bands  of  particles  may  be  the  result  of  processing 
an  MMC  where  the  original  material  consisted  of  clustered 
groups  of  particles  [Ref.  3]. 

Kocks  [Ref.  4]  showed  that  the  average  area  a  swept  out  by 
a  single  dislocation  segment  is  dependent  on  the  applied 
stress  a.  This  relationship  is  shown  in  Figure  1.  It  can  be 
seen  that  there  exist  a  stress  in  which  a  dislocation  can 
sweep  out  an  area  and  keep  moving  indefinitely.  This 
asymptotic  stress  is  the  macroscopic  flow  stress  a.  It  is 
also  apparent  that  the  greatest  flow  stress  for  a  given 
particle  density  is  for  a  regular  array  of  particles,  is  least 
for  a  clustered  array,  and  the  flow  stress  for  a  random  array 
lies  between  the  two.   Thus  the  distribution  of  second  phase 


particles  has  a  significant  role  in  determining  the  mechanical 
properties  of  two  phase  materials. 


a 

A 


Regular 


Clustered  (qualitative) 


I 

Figure  1.  Applied  stress  a  vs  mean  free  slip  area  a  of 
one  dislocation  loop,  for  strong  obstacles  of  area  density 
Ijl1,    in  square,  random  and  clustered  arrays  [Ref.  3]. 


Clearly,  it  would  be  desirable  to  have  a  means  of 
accurately  assessing  the  different  types  of  microstructures 
resulting  from  various  processing  techniques.  With  this 
knowledge  these  techniques  can  be  optimized  to  efficiently 
achieve  the  desired  properties  of  the  material.  Current 
methods  used  to  determine  these  values  assume  that  the 
distribution   is   random   and   that   the   particles   are 


infinitesimal  points,  which  is  certainly  not  the  case  for  real 
materials.  Methods  do  exist  for  finite  sized  particles  but 
here  the  assumption  is  that  they  are  distributed  in  a  regular 
geometric  array. 

This  study,  which  is  essentially  follows  on  from  that  of 
M.  Pas  [Ref.  5],  was  undertaken  to  guantify  such  values  as 
average  first  nearest  neighbor  particle  spacings  (5)  for 
various  area  fractions,  particle  size  distributions  and  for 
non-random  cases  as  well.  Any  trends  or  patterns  in  the 
associated  data  was  investigated  as  well.  In  order  to 
facilitate  mathematical  and  computer  analysis  computations 
were  kept  on  a  two  dimensional  basis. 


II.    BACKGROUND 

The  nearest  neighbor  separation  distance  between  second- 
phase  particles  governs  many  properties  of  two  phase 
materials.  Properties  such  as  strength,  ductility,  and 
toughness  are  sensitive  to  the  mean  separation  distance  and 
may  also  be  affected  by  the  range  and  uniformity  of  this 
distance  throughout  the  material.  Current  methods  for 
determining  this  spacing  for  particles  of  finite  size  is 
actually  based  on  the  assumption  that  the  particles  are 
infinitesimal  points  which  do  not  occupy  space.  Throughout 
this  thesis,  unless  otherwise  specified,  the  word  "point" 
shall  be  regarded  as  an  inf initesimally  small  entity  which 
does  not  occupy  space.  The  word  "particle"  shall  be  defined 
as  a  small  two  or  three  dimensional  object. 

Underwood  [Ref.  6]  has  described  a  method  for  determining 
the  average  first  nearest  neighbor  distance,  S,  for  randomly 
distributed  points.  The  method  is  based  on  the  probability  of 
finding  a  nearest  neighbor  to  a  given  point  within  an  annulus, 
of  radius  r  and  width  dr ,  centered  on  the  point  (Figure.  2) . 
This  probability  can  be  expressed  as: 


P(r)dr  = 


-   (rP(r) 
Jo 


dr 


2TzNadr  (2.1) 


'A 


where  NA  is  the  number  of  particles  per  unit  area  and  r  is  the 
distance  between  a  point  and  its  first  nearest  neighbor.  The 


Figure  2.  The  probability  of  locating  a  nearest  neighbor 
to  a  point  may  be  found  by  integrating  over  an  annulus 
centered  on  the  point. 


function  P(r)   which  satisfies  equation  (2.1)  was  shown  to  be 


P(r)    =  2nrNAe~*z2"* 


(2.2) 


Substituting  equation  (2.2)  into  the  expression  that  defines 
the  average  first  nearest  neighbor  distance  results  in  the 
following  expression  for  6  for  a  random  distribution  of 
infinitesimal  points: 


6  =  f~zp(z)di  =  o.sootf;1 

Jo 


-1/2 


(2.3) 


On  the  other  hand,  if  particles  or  points  are  arranged  in 
a  regular  hexagonal  array,  Figure  3,  the  nearest  neighbor 
distance  for  each  particle  will  be  the  same  everywhere  and 
will  be  equal  to  <SHex  for  the  array.  A  geometrical  analysis 
results  in  the  following  expression  for  5Hex  : 


6Hex.  =  1.07  5AT; 


-1/2 


(2.4) 


Figure  3.   Particles  arranged  in  a  hexagonal  array. 


Based  on  equations  (2.3)  and  (2.4)  it  can  be  seen  that, 
for  a  given  value  of  NA,  SHex  is  greater  than  S  for  a  random 
distribution. 

If  particles  (or  points)  are  arranged  into  any  regular 
array  the  first  nearest  neighbor  distance  (NND) ,  r,  would  be 
the  same  for  each  particle  and  therefore  would  equal  the 


resulting  value  of  S  for  the  array.  If  any  one  of  these 
particles  (or  points)  were  shifted  out  of  its  position  in  the 
array  it  would  come  closer  to  at  least  one  other  particle. 
Now,  there  would  be  at  least  two  particles  with  a  lower  first 
NND,  but  none  with  a  larger  first  NND  since  the  remaining 
particles  would  still  have  other  neighbors  at  a  distance  of  r. 
Therefore,  as  the  particles  are  shifted  out  of  a  perfect 
geometric  array  and  become  more  random  the  value  of  S  will 
decrease  towards  a  lower  limit  defined  by  the  Poisson 
distribution. 

Distributions  of  finite  size  particles  can  not  truly  be 
random.  The  area  (or  volume)  occupied  by  a  particle  removes 
the  possibility  that  another  particle  can  be  positioned  in  a 
portion  of  that  area.  This  is  the  effect  of  non-overlapping 
particles.  However,  for  dilute  arrays  (low  area  fractions) 
the  possibility  of  two  particles  overlapping  is  small. 
Therefore,  it  is  reasonable  to  expect  that  the  value  of  S  for 
dilute  arrays  of  finite  size  particles  would  in  fact  approach 
that  of  the  Poisson  distribution. 

In  previous  work,  performed  by  Mike  Pas  [Ref.  5],  attempts 
to  demonstrate  that  dilute  arrays  approached  the  Poisson 
distribution  were  not  completely  successful.  Pas  plotted  NA 
vs  S  for  various  area  fractions  of  constant  size  particles. 
The  plots  asymptotically  approached  values  of  S  greater  than 
those  of  the  Poisson  distribution  with  decreasing  area 


fraction.  However,  Pas  did  show  that  as  the  area  fraction  of 
particles  increases,  values  for  6  shift  towards  values 
approaching  that  of  a  hexagonal  array,  equation  (2.4).  Pas 
also  found  that  in  non-random  (high  and  low  density  bands  of 
particles)  distributions,  as  the  percentage  of  particles  in 
the  high  density  bands  increases,  the  overall  value  for  6 
decreases. 

The  ultimate  goal  of  this  work  would  be  an  analytical 
solution  to  quantify  the  degree  of  randomness  of  the  particle 
distributions  in  microstructures.  Incorporating  this  solution 
with  a  computer  based  image  analyzing  system  would  make  it 
possible  to  readily  assess  the  mechanical  properties  dependent 
on  the  particle  distributions. 
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III.     EXPERIMENTAL  PROCEDURE 

A.  OVERVIEW 

The  simulations  of  particle  distributions  were  generated 
using  a  computer.  Four  independent  programs  were  written,  one 
for  each  of  four  cases  analyzed:  random  distribution  of 
constant-size  particles  (RANPART1)  ;  non-random  distribution  of 
constant  size  particles  (RANPART2) ;  random  distribution  of 
particles  with  a  log-normal  size  distribution  (RANPART3) ;  and 
non-random  distribution  of  particles  with  a  log  normal  size 
distribution  (RANPART4) .  RANPARTl,  2,  3  and  4  are  included  in 
appendices  A,  B,  C  and  D  respectively. 

For  the  first  two  cases  the  algorithms  were  based  on 
programs  originally  developed  by  M.  Pas  [Ref.  5]  and  written 
in  TURBO  PASCAL.  Here,  Pas's  programs  were  rewritten  in 
FORTRAN  77  and  incorporated  several  other  modifications  as 
deemed  necessary.  The  programs  were  run  either  on  an  IBM- 
compatible  PC  or  a  Digital  Eguipment  Corp.  VAX  3100  work 
station. 

B.  RANDOM  DISTRIBUTION  OF  CONSTANT  SIZE  PARTICLES,  RANPARTl 

The  flow  chart  for  RANPARTl  is  shown  in  Figure  4.  In  the 
input  portion  of  the  program  the  user  enters  the  following 
parameters: 
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NO 


INPUT  INITIAL  PARAMETERS 


GENERATE  PARTICLE 
POSITION 


PARTICLE  OVERLAPS 
OTHER  PARTICLES 


NO 


STORE  PARTICLE  POSITION 


ALL  PARTICLES  POSITIONED 


YES 


SEGREGATE  PARTICLES 
ALONG  THE  BORDER 


DETERMINE  NEAREST 

NEIGHBOR  SPACING 

VALUES 


OUTPUT  DATA 


Figure  4.    Flow  chart  for  the  program  RANPARTlJOR 
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•  Area  percentage  of  particles  for  the  defined  region; 

•  x  dimension  of  the  defined  region; 

•  y  dimension  of  the  defined  region; 

•  Minimum  spacing  between  particles   (if  particles  are 
allowed  to  come  in  contact,  enter  0) ; 

•  Enter  "n"  if  a  specified  number  of  particles  is  required 
or  enter  "r"  if  a  specified  particles  radius  is  required; 

•  If  entered  "n" ,    enter  the  total  number  of  particles,  if 
entered   "r",  enter  the  particles  radius; 

•  Enter  "A"  if  particles  are  of  finite  size  or  enter  "p"  if 
the  particles  are  infinitesimal  points. 

From  these  parameters  the  total  area  of  the  defined  region, 

total  area  occupied  by  the  particles,  the  area  per  particle, 

and  the   total   number   of   particles   or  particle   radius 

(whichever  value  was  not  input)  are  calculated. 

The  particles  are  positioned  in  the  next  portion  of  the 
program.  A  pseudorandom  number  generator  is  called  twice  to 
provide  x  and  y  coordinates  for  a  given  particle.  This 
particle  is  then  checked  to  ensure  that  it  does  not  overlap 
any  other  previously  sited  particle.  This  is  accomplished  by 
calculating  the  distance  to  every  other  particle.  If  this 
distance  is  less  than  twice  the  sum  of  the  particle  radius  and 
the  minimum  spacing,  the  position  is  rejected  and  the  program 
loops  back  to  generate  a  new  set  of  coordinates.  This  process 
is  repeated  until  the  total  number  of  particles  required  to 
establish  the  desired  area  fraction  has  been  met. 

Once  all  the  particle  locations  have  been  established, 
particles  located  along  the  border  are  segregated  from  those 
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in  the  interior  of  the  defined  region.  This  is  done  so  that 
errors  in  first  nearest  neighbor  distance  (NND)  calculations 
associated  with  the  border  will  be  minimized.  The  underlying 
rationale  for  this  step  in  the  program  will  be  discussed  in 
the  next  chapter.  All  particle  positions  are  checked  and 
those  having  coordinates  placing  them  within  one  particle 
radius  of  the  border  are  moved  to  a  separate  array  in  the 
program  memory. 

First  NND  calculations  are  made  by  calculating  for  each 
interior  particle  the  distance  to  every  other  particle, 
including  those  located  along  the  border,  and  then  finding  and 
saving  the  minimum  first  NND  value.  This  process  is  repeated 
for  every  other  interior  particle.  The  average  first  NND,  S, 
is  calculated  by  summing  the  minimum  first  NND  values  and 
dividing  by  total  number  of  interior  particles  . 

In  order  to  plot  the  distributions  of  first  NND  values  an 
increment  had  to  be  determined.  To  do  this,  the  Sturgis  rule 
was  applied  as  it  is  an  established  method  for  determining  the 
number  of  class  intervals  from  a  set  of  values.  The  Sturgis 
rule  is  defined  by  the  following  expression: 

N  =   1  +  3.31og(fl)  (3.1) 

Where  N  is  the  number  of  class  intervals  and  n  is  the 
number  of  values  in  the  set.  The  value  of  n  for  this  work  is 
equal  to  the  number  of  interior  particles.  The  number  of 
class  intervals  was  divided  into  the  difference  between  the 


14 


maximum  and  minimum  first  NND  values.  This  result  (defined  as 
"dpinc"  in  the  programs)  was  used  as  the  incremental  range  for 
sorting  first  NND  values  into  equal  sized  intervals.  The 
number  of  NND  values  in  each  interval  was  then  divided  by  the 
total  number  of  interior  particles  to  get  a  relative  quantity. 
The  output  of  the  program  consists  of  all  input  parameters, 
the  calculated  values  described  above  and  a  calculated  value 
for  NA,  the  number  of  particles  per  unit  area. 

RANPART1  was  run  for  area  fractions  ranging  from  1%  to  50% 
utilizing  a  particle  radius  of  0.5  units.  The  defined  area 
was  adjusted  in  each  run  to  maintain  a  nominal  value  Np~5000 
particles.  The  minimum  spacing  was  kept  at  0  for  all 
simulations. 

Each  run  was  repeated  using  the  same  parameters,  but  now 
declaring  the  particles  to  be  points.  The  results  consisted 
of  point  distribution  and  a  distribution  of  finite  sized 
particles  with  the  same  value  of  NA. 

For  each  simulation  equation  (2.1)  was  used  to  plot  a 
Poisson  distribution  corresponding  the  simulated  distribution 
of  particles  and  points.  The  value  of  NA  came  from  the 
simulated  point  distribution.  The  probability  density 
function  was  plotted  as  a  function  of  r  in  the  range  from  the 
minimum  and  maximum  values  of  first  NND  derived  from  the  point 
distribution.  The  calculated  increment  Ar,  based  on  the 
Sturgis  rule,  was  substituted  into  the  equation  for  dr. 
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C.   NON-RANDOM   DISTRIBUTION   OF   CONSTANT   SIZE   PARTICLES, 
RANPART2 

The  following  modifications  were  made  to  RANPART1  to 
produce  RAN PART 2 .  The  defined  area  in  which  the  particles  are 
positioned  was  divided  into  five  sub-areas  of  equal  size.  In 
the  input  portion  of  the  program  the  user  declares  what 
proportion  of  the  total  area  occupied  by  the  particles  will 
be  assigned  to  each  subsection.  The  sum  of  the  percentages 
must  equal  100%. 

After  each  particle  (for  finite  sized  particles)  is 
checked  for  overlap,  the  program  determines  which  sub-area  the 
particle  (or  point)  belongs.  By  maintaining  a  running  sum  on 
the  area  occupied  by  the  particles  in  each  sub-area,  the 
computer  determines  whether  the  addition  of  another  particle 
in  a  given  sub-area  would  cause  the  area  occupied  by  the 
particles  to  exceed  the  limit.  The  particle  position  is 
rejected  if  the  limit  will  be  exceeded.  If  this  happens  the 
program  loops  back  to  generate  another  pair  of  particle 
coordinates.  When  the  process  is  complete  the  desired  degree 
of  banding  has  been  achieved.  Figure  5  is  the  flow  chart  for 
RANPART2 . 

All  runs  made  for  RANPART2  used  10%,  35%,  10%,  35%,  and 
10%  sub-area  particle  densities.  This  was  considered  to 
provide  a  relatively  large  deviation  between  the  high  and  low 
density  bands.  The  simulations  were  made  for  area  fractions 
ranging  from  0.01  to  0.30. 
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Figure5.  Flow  chart  for  the  program  RANPART2.FOR 
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D.   RANDOM  DISTRIBUTION  OF  PARTICLES  WITH  A  LOGNORMAL  SIZE 
DISTRIBUTION,  RANPART3 

In  RANPART3  a  second  random  number  generator  is  used  to 
determine  a  particle  radius  from  a  lognormal  distribution. 
The  call  for  the  particle  radius  is  made  before  the  call  for 
the  particle  coordinates  so  that,  if  particle  overlap  occurs, 
only  a  new  position  will  be  generated  but  the  existing  radius 
will  still  be  used.  Only  when  a  particle  is  finally 
positioned  or,  if  particle  overlap  occurs  600  times 
consecutively  for  the  same  particle  radius,  is  a  new  particle 
radius  generated.  This  process  is  required  in  order  to 
maintain  a  lognormal  size  distribution  and  to  prevent  biasing 
towards  smaller  sized  particles  as  the  particle  density 
increases.  The  program  must  now  use  individual  particle  radii 
in  subsequent  calculations. 

The  input  portion  of  the  program  was  modified  as  follows. 
The  user  must  enter  a  value  approximately  equal  to  the  natural 
logarithm  of  the  mean  particle  radius  desired,  which  is  the 
mean  of  the  underlying  normal  distribution,  and  a  standard 
deviation.  Several  iterations  of  these  two  variables  may  be 
required  to  get  the  desired  mean  radius,  especially  for  larger 
standard  deviations.  Since  a  point  distribution  derived  from 
the  same  input  parameters  of  constant  size  particles  would  be 
no  different  than  that  for  particles  of  varying  size,  this 
option  is  not  available  for  the  programs  with  a  lognormal 
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particle  size  distribution.  The  total  number  of  particles 
entered  in  the  input  portion  of  the  program  must  be  of 
sufficient  quantity  to  ensure  that  the  number  of  interior 
particles  generated  is  enough  to  meet  the  desired  area 
fraction.)  Since  now  the  particles  are  no  longer  of  constant 
size  there  is  no  direct  relationship  between  the  number  of 
particles  and  the  total  area  occupied  by  the  particles.  A 
running  sum  of  calculated  particle  areas  is  maintained. 
Particle  generation  continues  until  this  sum  is  equal  to  the 
value  corresponding  to  the  desired  area  fraction. 

Segregation  of  particles  along  the  border  is  performed 
much  the  same;  but  here  the  largest  particle  radius  in  the 
distribution  is  used  to  establish  the  extent  of  the  border. 
Any  particle  whose  center  lies  outside  this  border  is 
segregated. 

The  maximum,  minimum  and  mean  particle  radius  are 
determined  and  provided  at  the  output.  The  flow  chart  for 
RANPART3  is  shown  in  Figure  6. 

RANPART3  was  run  for  the  same  area  fractions  as  RANPART1. 
Values  of  -3.07  and  0.4  were  used  as  input  for  the  mean  of  the 
underlying  normal  distribution  and  for  the  standard  deviation, 
respectively. 
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Figure  6.  Flow  chart  for  the  program  RANPART3.FOR 
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E.   NON-RANDOM  DISTRIBUTION  OF  PARTICLES  WITH  A  LOGNORMAL  SIZE 
DISTRIBUTION,  RANPART4 

RANPART4  was  derived  from  RANPART3  in  the  same  way 
RANPART2  was  derived  from  RANPART1.  Five  sub-areas  were 
incorporated  into  RANPART4  just  as  they  were  incorporated  into 
RANPART2 .  To  maintain  a  record  of  the  area  occupied  by  the 
particles  in  each  sub-area,  the  area  of  each  new  particle 
needs  to  be  calculated  and  added  to  the  existing  area  in  which 
it  lies. 

The  flow  chart  for  RANPART4  is  shown  in  Figure  7. 
RANPART4  was  run  for  the  same  area  fractions  as  RANPART2  and 
for  the  same  values  of  mean  and  standard  deviation  as 
RANPART3 . 
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Figure  7.  Flow  chart  for  the  program  RANPART4.FOR 
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IV.    RESULTS  AND  DISCUSSION 

A.  OVERVIEW 

The  first  phase  of  this  work  consisted  of  a  study  to 
determine  a  statistically  sufficient  number  of  particles  to  be 
used  in  the  simulations.  The  problem  of  particles  positioned 
along  the  border  of  the  array  was  also  addressed  in  this 
phase.  In  the  second  phase,  four  general  cases  were  studied: 
random  distribution  of  points  and  finite  sized  particles  of 
constant  size;  non-random  distribution  of  points  and 
particles  of  constant  size;  random  distribution  of  particles 
with  a  log  normal  size  distribution;  and  finally,  non-random 
distribution  of  particles  with  a  log  normal  size  distribution. 

B.  DETERMINATION  OF  A  SUFFICIENT  NUMBER  OF  PARTICLES  AND  THE 
PROBLEM  OF  PARTICLES  POSITIONED  ALONG  THE  BORDER 

In  order  to  determine  a  statistically  sufficient  number  of 
particles  to  be  used  in  the  simulations  the  program  RANPART1 
was  modified.  The  modifications  essentially  consisted  of 
enclosing  RANPART1  inside  two  loops.  This  allowed  repeated 
program  running  for  the  same  parameters.  Then,  the  average  of 
the  results  of  repeated  runs,  was  calculated.  This  experiment 
was  used  to  calculate  an  average  value  of  S,  over  a  given 
number  of  runs,  for  various  numbers  of  particles  (Np)  , 
holding  NA  constant.   The  test  was  conducted  for  two  area 
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fractions,  0.10  and  0.40  and  for  point  distributions  as  well. 
NA  was  arbitrarily  set  equal  to  100.  The  results  of  the  tests 
are  shown  in  Table  1  and  in  Figure  8. 


Np 

af  =  0.1 

af  =  0.4 

Points 

No. 
Runs 

10 

0.06801 

0.08902 

0.05781 

1000 

20 

0.06562 

0.08633 

0.05556 

1000 

30 

0.06442 

0.08524 

0.05425 

1000 

40 

0.06376 

0.08455 

0.05362 

1000 

50 

0.06349 

0.08434 

0.05326 

1000 

60 

0.06318 

0.08400 

0.05282 

500 

70 

0.06282 

0.08370 

0.05271 

500 

80 

0.06281 

0.08352 

0.05257 

500 

90 

0.06241 

0.08348 

0.05240 

500 

100 

0.06233 

0.08326 

0.05249 

100 

200 

0.06173 

0.08255 

0.05146 

100 

300 

0.06154 

0.08223 

0.05131 

100 

400 

0.06120 

0.08211 

0.05106 

100 

500 

0.06127 

0.08198 

0.05107 

50 

600 

0.06107 

0.08195 

0.05083 

50 

700 

0.06099 

0.08181 

0.05098 

25 

800 

0.06091 

0.08172 

0.05046 

25 

900 

0.06087 

0.08172 

0.05087 

25 

1000 

0.06076 

0.08178 

0.05075 

5 

2000 

0.06059 

0.08150 

0.05025 

5 

3000 

0.06035 

0.08130 

0.05059 

5 

4000 

0.06048 

0.08130 

0.05035 

5 

5000 

0.06041 

0.08127 

0.05030 

5 

Table  1.   Results  from  the  modified  RANPART1  for  area 
fractions  of  0.1  and  0.4  and  for  points,  NA  =  100. 
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Figure    8.      Total    number    of    particles    vs    average    first    nearest 
neighbor    distance    for    Na    "    100 


25 


For  the  random  distribution  of  points,  S  approaches  the 
theoretical  value  of  0.05  for  NP  >  1000.  For  both  area 
fractions,  0.40  and  0.10,  S  approaches  and  appears  to  remain 
constant  at  values  of  about  0.081  and  0.060  respectively, 
again  at  Np  >  1000  in  both  cases. 

In  developing  the  method  for  positioning  circles  or  points 
on  a  plane  the  dimensions  of  the  planar  surface  had  to  be 
established.  Within  this  finite  sized  area  all  the  particles 
were  positioned.  From  this  situation  it  became  evident  that 
when  calculating  first  nearest  neighbor  distances,  a  particle 
whose  center  lies  on  or  near  the  border  of  the  defined  region 
could  conceivably  have  a  closer  neighbor  than  that  calculated. 
Figure  9  shows  a  defined  area  and  a  particle,  a,  located  near 
the  border.  Particle  a's  nearest  neighbor  within  the  area  is 
particle  b.  However,  it  is  possible  that  if  the  defined  area 
had  been  larger,  particle  a  may  have  a  nearest  neighbor  even 
closer  than  b  such  as  particle  c.  Clearly  this  would  be  the 
situation  in  any  real  material  (imagine  the  field  of  view  of 
a  micrograph  being  expanded) .  Since  the  existence  of  these 
additional  particles  could  only  serve  to  lower  the  first  NND, 
any  calculations  made  for  S  without  taking  this  into 
considerations  would  be  overestimated. 

Calculations  for  6  were  repeated,  holding  NA  constant  and 
for  the  same  area  fractions,  as  described  above.  The  model 
was  modified  such  that  first  NND  calculations  were  not  made 
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Figure  9.  The  border  problem.  Particles  positioned  along  the  border  of 
a  defined  area  (a)  may  have  a  potential  nearest  neighbor  (c),  if  the 
defined  region  was  expanded,  which  is  closer  than  the  nearest  neighbor 
positioned  inside  the  area  (b). 
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for  particles  whose  centers  fell  within  one  particle  radius 
from  the  border.  These  border  particles  were  used  in 
determining  first  NND  for  particles  located  inside  this 
exclusion  range.   The  results  are  shown  in  Figure  10. 

This  attempt  to  eliminate  the  border  problem  was  only 
partially  successful  as  shown  in  Figure  10  and  has  the 
greatest  effect  for  low  numbers  of  particles.  The  values  for 
6  were  consistently  lower  for  calculations  involving  both 
particles  and  points  when  the  numbers  of  points  or  particles 
utilized  was  below  1000.  For  NP  greater  than  1000  the  effects 
of  using  a  larger  number  of  particles  or  points  compensates 
for  the  border  problem  by  simply  increasing  the  relative 
number  of  particles  located  inside  the  defined  area. 

The  data  compiled  in  Table  1  was  used  to  determine  a  more 
exact  solution  for  the  asymptotic  value  of  S  (as  Np  -»  oo)  and 
show  that  it  is  a  function  of  the  5-^^  and  Np.  Close 
examination  of  Figure  8  reveals  the  similar  form  of  all  three 
curves.  That  is,  for  any  value  of  Np  the  difference  between 
<W  and  *Theo.  (for  points)  or  <Sobs  and  <SAsym  (for  particles) 
appears  to  be  the  same.  Based  on  the  previous  observation 
that  this  difference  is  due  to  the  particles  positioned  along 
the  border,  the  following  relationship  was  defined: 

^A6  =6ois.  -  6rheo.        (4.1) 
where  NpB  is  the  number  of  particles  along  the  border  and  A<5 
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Figure    10.       Total    number    of    particles    vs    average    first    nearest 
neighbor    distance    for    Na    "    100.      The    three    pairs    of    curves 
demonstrate    the    affect    of    particles    located    on    the    border. 
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is  the  overestimate  in  S  associated  with  each  particle  along 
the  border.  For  a  random  distribution  of  particles  in  a 
square  region,  Np  can  be  taken  to  be  equal  to  the  square  of 
the  number  of  particles  along  the  border,  Nj,,2.  Therefore  N^ 
=  Np1/2  and  equation  (4.1)  can  be  expressed  as: 

AS  =  bobs.   -  bTheo.  (4.2) 


np 

Figure  11  is  a  plot  of  Np"1/2  vs  <Sot)S  -  5^^  from  the  data 
compiled  in  Table  1.  These  curves  clearly  define  a  linear 
relationship  from  which  the  following  equations  can  be 
derived: 

Points:      60bs_    =  0.0235-^    +  6Theo_  (4.3) 


af  =  0.1:      60jbs     =  0.0248  —    +  0.0100    +  6rheo  (4.4) 


af=0.4:      bobs_    =  0.0248-^    +  0.0309    +  6Theo.  (4.5) 

By  taking  the  limit  as  Np  -►  <»  in  each  equation  (4.3),  (4.4) 
and  (4.5),  <50bs  -♦  5Asym..  The  solutions  for  <SAsym  in  these 
equations  are  0.0500,  0.0600  and  0.0809,  respectively.  In 
general,  this  value  for  <SAsym  can  be  taken  to  be  the  average 
first  nearest  neighbor  distance  for  any  distribution  of  points 
or  particles  of  finite  size. 


30 


0.04 


0.03 


o 

H 

1      0.02      — 


CD 

O 
«o 


0.01 


0.00 


6Theo.  =    0.0309    +    0.0248Np 


a  Particles,  a*  =  0.40,  rp  =  0.05 
a  Particles,  af  =  0.10,  rp  =  0.05 
o    Calculated   Points 


Theo.    =    0.0100    +    0.0248Np 


0.5 


Obs.  "     °Theo.    =    0.0234NP 


-0.5 


0.00  0.05  0.10  0.15  0.20  0.25  0.30  0.35 


Np 


0.5 


-0.5 


&Ure   J1*  ,?P       vs    6Obs.     "     5Theo.for    area    factions    of 
10    and    0.40. 


31 


Based  on  the  results  above  all  subsequent  simulations  were 
made  using  both  the  border  particle  exclusion  algorithm  and  a 
nominal  value  of  5000  particles  or  points. 

C.   RANDOM  DISTRIBUTION  OF  CONSTANT  SIZE  PARTICLES 

For  all  values  of  NA  evaluated,  1.275  <  NA  <  63.82,  the 
calculated  distribution  for  random  points  very  closely 
approximated  the  corresponding  Poisson  distribution  (Figure  12 
-16).  For  area  fractions  less  than  0.05,  as  shown  in  Figure 
12,  the  particle  distributions  generated  approximated  the 
Poisson  as  well.  In  these  dilute  distributions  the  particles 
behave  like  points  since  the  likelihood  of  two  randomly 
generated  particles  overlapping  is  not  significant.  The 
difference  between  the  distributions  is  that  the  particle 
distribution  becomes  truncated  at  a  value  equal  to  twice  the 
particle  radius.  This  was  a  result  of  the  non-overlap 
criterion  preventing  the  particle  centers  from  getting  closer 
than  one  particle  diameter.  Figures  17  -  21  are  the 
corresponding  plots  of  the  actual  particle  distributions 
represented  in  Figures  12  -  16. 

Figure  2  2  demonstrates  the  effect  that,  as  the  area 
fraction  is  increased,  the  distribution  is  shifted  towards 
lower  values  of  NND.  A  shift  in  S  from  the  theoretical 
Poisson  value  (5-rheo.  =  0.5NA1/2)  towards  the  value  for  a 
hexagonal  array  (6Hex  =  1.075NA"1/2)  was  also  observed.  At  high 
area  fractions  the  distributions  for  particles  approach  delta 
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PARTICLE  AREA  FRACTION:  0.01 

PARTICLE  RADIUS:  0.05 

AVERAGE  NEAREST  NEIGHBOR  SPACING:  0.4608 

NUMBER  PARTICLES  PER  UNIT  AREA:  1.273 

Figure  17.   Random  distribution  of  constant  size  particles. 
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Figure  18.    Random  distribution  of  constant  size  particles. 
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Figure  19.    Random  distribution  of  constant  size  particles. 
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Figure  20.    Random  distribution  of  constant  size  particles. 
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AVERAGE  NEAREST  NEIGHBOR  SPACING:  0.1082 

NUMBER  PARTICLES  PER  UNIT  AREA:  62.82 


Figure  21.   Random  distribution  of  constant  size  particles. 
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functions  with  the  minimum  first  NND  having  the  highest 
frequency  of  occurrence. 

The  effects  of  increasing  area  fraction  for  finite  sized 
particles  is  shown  in  Figure  23.  As  NA  increases  from  values 
of  very  diluted  arrays  to  more  dense  distributions  the 
corresponding  values  of  5  shift  from  that  of  a  random 
distribution  of  infinitesimal  points  to  that  of  a  hexagonal 
array.  Figure  2  3  shows  that  this  trend  is  consistent  for 
various  particle  radii.  The  reader  should  note  the  asymptotic 
approach  towards  the  random  theoretical  distribution  of 
points.  This  is  indicative  of  using  a  sufficient  number  of 
particles  in  the  simulations. 

D.   NON-RANDOM  DISTRIBUTION  OF  CONSTANT  SIZE  PARTICLES 

A  non-random  distribution  was  simulated  by  dividing  the 
total  defined  area  in  to  five  sub-areas  of  equal  size.  In  all 
simulations  there  were  three  low  density  zones  separated  by 
two  high  density  zones.  The  high  density  zones  were  each 
assigned  3  5%  of  the  particles  occupying  the  total  area,  while 
the  low  density  zones  were  each  assigned  10%.  These 
parameters  were  used  for  all  simulations. 

The  non-random  distribution  of  constant  size  particles  for 
various  area  fractions  are  shown  in  Figures  24  -  27.  Plots  of 
the  particles  themselves  are  shown  in  Figures  28  -  31.  These 
figures  were  compared  with  the  corresponding  distributions  and 
particle  plots,  with  respect  to  area  fraction,  in  the  random 


44 


10 


m 
U 
2 
< 

H 
V3 


o 

so 
K 
O 

M 

w 
2 

H 

a: 
< 


H 


W 

U 

< 

w 

> 


1     — 


.01  .1  1  10  100 

NUMBER    OF    PARTICLES    PER    UNIT    AREA 

Figure    23.        Number     of    particles    per     unit    area    vs    average    first 
nearest    neighbor    distance    for    particles    of    varying    radius. 
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Figure  28.    Non— random  distribution  of  constant  size  particles. 
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Figure  29.    Non-random  distribution  of  constant  size  particles. 
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Figure  30.    Non— random  distribution  of  constant  size  particles. 
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Figure  31.    Non-random  distribution  of  constant  size  particles. 
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distribution  case.  It  was  observed  that  the  fundamental  form 
of  the  distributions,  trend  in  values  of  the  observed  S,  and 
first  NND  were  similar. 

Additionally,  in  the  non-random  case,  there  was  a  small 
but  consistent  shift  in  the  distribution  towards  lower  values 
of  first  NND.  This  pattern  was  attributed  to  the  presence  of 
the  high  density  bands.  For  a  given  NA,  the  greater  degree  of 
packing  of  points  or  particles  in  the  high  density  zones 
caused  the  first  NND  to  shift  toward  lower  values.  This 
effect  was  also  reflected  in  a  greater  frequency  of  the 
minimum  first  NND  when  compared  to  that  of  the  random  case. 
A  summary  of  the  observed  trends  in  the  non-random/constant 
particle  size  distributions  are  shown  in  Figure  32. 

The  corresponding  particle  plots,  comparing  the  relative 
densities  of  the  high  and  low  density  bands,  are  shown  in 
Figures  28  -  31.  Compare  these  plots  with  those  of  the  random 
distributions  for  the  same  area  fractions  (Figures  17  -  21) . 
For  af  =  0.30  the  difference  between  the  random  (Figure  20)  and 
non-random  (Figure  31)  distributions  is  completely  obvious  to 
the  eye.  There  is  a  clear  distinction  between  the  high  and 
low  density  bands  in  the  non-random  distribution.  This 
distinction  makes  it  easy  for  the  observer  to  tell  the 
difference  between  the  random  and  non-random  particle 
distributions  based  on  these  plots.  Again,  for  af  =  0.30, 
compare  the  distribution  curves  for  the  random  (Figure.  15) 
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and  non-random  (Figure  27)  cases.  The  distinction  is  not 
great  and  the  discernable  difference  is  a  greater  maximum 
frequency  for  the  lowest  first  NND  in  the  non-random  case. 

The  visual  difference  in  the  particle  pl6ts  is  much  more 
subtle  for  low  area  fractions  compared  to  high  area  fractions. 
Compare  the  random  and  non-random  particle  plots  for  af  =  0.01 
(Figures  17  and  28)  .  Unlike  the  plots  for  the  30%  area 
fraction,  the  distinction  is  not  obvious.  However,  the  same 
can  be  said  for  the  NND  distribution  curves  (Figures  12  and 
24)  .  For  the  low  area  fraction  both  curves  approach  the 
distribution  of  random  points.  This  made  any  distinction 
between  them  even  more  subtle  than  the  curves  for  the  3  0%  area 
fraction. 

E.   RANDOM  DISTRIBUTION  OF  PARTICLES  WITH  A  LOGNORMAL  SIZE 

DISTRIBUTION 

In  actual  microstructures  the  second  phase-particle  size 
generally  will  not  be  constant  throughout  the  material.  In 
fact  it  has  often  been  reported  that  second-phase  particles 
follow  a  lognormal  size  distribution.  Such  was  the  basis  for 
the  following  experiments. 

The  same  area  fractions  were  analyzed  for  particles  with 
a  lognormal  size  distribution  (Figures  33  -  37)  as  with  those 
for  constant  size  particles.  The  corresponding  particle  plots 
are  shown  in  Figures  38  -  42.  For  very  dilute  arrays  of 
particles,  the  distributions  approached  very  closely  to  those 
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Figure   33.      Random    distribution   of   points   and   particles   with   a 
lognormal    size    distribution. 
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Figure    34.      Random    distribution    of    points    and    particles    with    a 
lognormal    size    distribution. 
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Figure   35.      Random    distribution   of    points    and    particles    with   a 
lognormal    size    distribution. 
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Figure    36.      Random    distribution    of    points    and    particles    with    a 
lognormal    size    distribution. 
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Figure   37.      Random    distribution    of    points   and    particles   with    a 
lognormal    size    distribution. 
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PARTICLE  AREA  FRACTION:  0.01 

MEAN  PARTICLE  RADIUS:  0.0499 

AVERAGE  NEAREST  NEIGHBOR  SPACING:  0.4944 

NUMBER  PARTICLES  PER  UNIT  AREA:  1.094 

Figure  38.    Random  distribution  of  particles  with  a  lognormal 
size  distribution. 
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Figure  39.    Random  distribution  of  particles  with  a  lognormal 
size  distribution. 
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PARTICLE  AREA  FRACTION:  0.10 

MEAN  PARTICLE  RADIUS:  0.0503 

AVERAGE  NEAREST  NEIGHBOR  SPACING:  0.1808 

NUMBER  PARTICLES  PER  UNIT  AREA:  10.75 

Figure  40.    Random  distribution  of  particles  with  a  lognormal 
size  distribution. 
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PARTICLE  AREA  FRACTION:  0.30 

MEAN  PARTICLE  RADIUS:  0.0506 

AVERAGE  NEAREST  NEIGHBOR  SPACING:  0.1256 

NUMBER  PARTICLES  PER  UNIT  AREA:  31.56 

Figure  41.    Random  distribution  of  particles  with  a  lognormal 
size  distribution. 
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PARTICLE  AREA  FRACTION:  0.50 

MEAN  PARTICLE  RADIUS:  0.0479 

AVERAGE  NEAREST  NEIGHBOR  SPACING:  0.1012 

NUMBER  PARTICLES  PER  UNIT  AREA:  59.65 

Figure  42.    Random  distribution  of  particles  with  a  lognormal 
size  distribution. 
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calculated  for  points  and  to  the  Poisson  distribution  (Figure 
33)  .  But  now,  as  the  area  fraction  is  increased  the  particle 
distribution  still  maintains  a  shape  similar  to  the  Poisson 
distribution.  Unlike  the  case  for  constant  size  particles, 
the  minimum  first  NND  can  be  smaller  than  the  mean  particle 
diameter.  A  smaller  particle  can  always  be  fitted  between  two 
or  more  larger  ones,  thus  reducing  the  first  NND  below  the 
value  of  the  mean  particle  diameter.  This  also  resulted  in  a 
more  nearly  constant  peak  frequency  of  first  NND  for  all 
ranges  of  area  fractions  examined. 

The  trends  for  6  of  the  particles  and  ranges  of  first  NND 
(Figure  43)  were  the  same  as  those  in  the  constant  size 
particle  case  for  lower  values  of  NA.  At  larger  values  of  NA 
the  value  of  6  for  lognormally  distributed  particles  becomes 
less  than  that  for  the  corresponding  constant  particle  sizes. 

F.   NON-RANDOM  DISTRIBUTION  OF  PARTICLES  WITH  A  LOGNORMAL  SIZE 

DISTRIBUTION 

The  same  parameters  for  banding  used  in  the  non-random, 
constant  size  particle  case  were  used  in  the  case  of  non- 
random  distribution  of  particles  with  a  lognormal  size 
distribution.  The  results  of  this  analysis  are  shown  in 
Figures  44  -  47  for  the  distributions  and  Figures  48  -  51  for 
the  particle  plots.  The  features  are  similar  in  many  respects 
to  those  reported  earlier  for  non-random/constant  particle 
size  distributions.   As  the  area  fraction  increases  S   shifts 
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Figure    44.      Non-random    distribution    of    points    and    particles    with    a 
lognormal    size    distribution. 
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Figure    45.      Non-random    distribution    of    points    and    particles    with    a 
lognormal    size    distribution. 
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Figure    46.      Non-random    distribution    of    points    and    particles   with    a 
lognormal    size    distribution. 
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Figure    47.      Non-random    distribution    of    points   and   particles   with   a 
lognormal    size    distribution. 
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TOTAL  PARTICLE  AREA  FRACTION:  0.01 

SUB-AREA  PARTICLE  AREA  FRACTIONS:  0.10,  0.35, 

0.10,  0.35,  0.10 

MEAN  PARTICLE  RADIUS:  0.0502 

AVERAGE  NEAREST  NEIGHBOR  SPACING:  0.4-703 

NUMBER  PARTICLES  PER  UNIT  AREA:  1.074 

Figure  48.    Non-random  distribution  of  particles  with  a  log- 
normal  size  distribution. 
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TOTAL  PARTICLE  AREA  FRACTION:  0.05 

SUB-AREA  PARTICLE  AREA  FRACTIONS:  0.10,  0.35, 

0.10,  0.35,  0.10 

MEAN  PARTICLE  RADIUS:  0.0503 

AVERAGE  NEAREST  NEIGHBOR  SPACING:  0.2290 

NUMBER  PARTICLES  PER  UNIT  AREA:  5.377 

Figure  49.    Non— random  distribution  of  particles  with  a  log- 
normal  size  distribution. 
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TOTAL  PARTICLE  AREA  FRACTION:  0.10 

SUB-AREA  PARTICLE  AREA  FRACTIONS:  0.10,  0.35, 

0.10,  0.35,  0.10 

MEAN  PARTICLE  RADIUS:  0.0508 

AVERAGE  NEAREST  NEIGHBOR  SPACING:  0.1762 

NUMBER  PARTICLES  PER  UNIT  AREA:  10.53 

Figure  50.    Non-random  distribution  of  particles  with  a  log- 
normal  size  distribution. 
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TOTAL  PARTICLE  AREA  FRACTION:  0.30 

SUB-AREA  PARTICLE  AREA  FRACTIONS:  0.10,  0.35, 

0.10,  0.35,  0.10 

MEAN  PARTICLE  RADIUS:  0.0470 

AVERAGE  NEAREST  NEIGHBOR  SPACING:  0.1100 

NUMBER  PARTICLES  PER  UNIT  AREA:  37.79 

Figure  51.    Non— random  distribution  of  particles  with  a  log- 
normal  size  distribution. 
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from  a  very  close  approximation  to  the  value  for  a  Poisson 
distribution  towards  that  of  a  hexagonal  array.  Unlike  the 
cases  for  constant  size  particles,  the  distribution  curve 
maintains  a  shape  similar  to  the  Poisson  distribution  and  does 
not  reach  a  peak  frequency  of  first  NND  corresponding  to  the 
minimum  first  NND.  This,  again,  is  due  to  the  possibility  of 
positioning  a  smaller  particle  between  two  or  more  larger 
ones.  There  is  a  consistent  increase  in  the  maximum 
occurrence  of  any  first  NND  value  over  that  of  the  random 
case.  This  is  also  a  result  of  the  greater  density  of 
particles  in  the  high  density  zones.  Figure  52  summarizes  the 
shift  in  the  distributions  as  the  area  fraction  increases. 

In  comparing  the  particle  plots  for  non-random 
distributions  of  particles  with  a  lognormal  size  distribution, 
Figures  4  8  -  51,  to  those  of  the  random  case,  the  same 
observations  are  made  as  for  the  constant  particle  size 
distributions.  For  relatively  high  area  fractions  the 
difference  between  random  and  non-random  particle  dispersions 
is  obvious  to  the  eye.  But  for  very  low  area  fractions  it 
becomes  more  difficult  to  differentiate  between  the  two.  This 
observation  is  summarized  graphically  in  Figures  53,  54  and 
55.  The  first  two  figures  represent  random  and  non-random 
distributions,  respectively,  for  particles  of  constant  and 
varying  size.  Both  figures  are  indeed  very  similar,  each 
showing  the  effect  of  smaller  particles  positioned  between 
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Figure    53.       Number    of    particles    per    unit    area    vs    average    first 
nearest    neighbor    distance;    for    random    distributions    of    constant 
size    particles    and    for    particles    with    a    lognormal    size    distribution, 
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Figure    54.      Number    of    particles    per    unit    area    vs    average    first 
nearest    neighbor    distance;    for    non-random    distributions    of    constant 
size    particles    and   for    particles   with    a    lognormal    size    distribution. 


80 


Random     Distribution 
Non-random    Distribution 


U 

< 

H 


O 

CO 

u 

w 

H 
(XI 

UJ 
OS 

< 

w 
2 

H 
CD 


(J 
< 

w 
> 


.1 


1  10  100 

NUMBER    OF    PARTICLES    PER    UNIT    AREA 

Figure    55.       Number    of    particles    per    unit    area    vs    average    first 
nearest    neighbor    distance;    for    random    and    non-random 
distributions    of    particles    with    a    log    normal    size    distribution, 
rm    -    0.05,    a   -   o.4. 
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larger  ones  reducing  the  value  of  S  for  a  given  NA.  Further, 
in  comparing  the  two  figures  the  minimal  deviation  between 
random  and  non-random  distributions  was  reflected  in  the  close 
similarities  of  the  curves.  This  observation  is  more  apparent 
in  figure  55  for  the  case  of  random  and  non-random 
distributions  of  particles  with  a  lognormal  size  variance. 

The  analysis  conducted  shows  that  the  human  eye  is  very 
capable  of  detecting  non-randomness  of  particle  distributions 
for  sufficiently  dense  arrays.  However,  when  presented  with 
the  statistical  data,  the  discernment  of  random  or  non-random 
is  not  so  obvious.  The  degree  of  non-randomness  imposed  on 
the  model  was  judged  to  be  relatively  high,  35%  and  10%  for 
high  and  low  density  bands  respectively.  Had  the  difference 
between  the  two  bands  been  less,  there  would  have  been  even 
greater  similarities  between  the  random  and  non-random 
distributions.  The  question  then  becomes;  is  it  possible  to 
develop  a  system  which  can  differentiate  between  random  and 
non-random  particle  distributions  and  quantify  the  results? 
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V.    CONCLUSIONS  AND  RECOMMENDATIONS 

The  following  conclusions  may  be  drawn  from  the  work 
described  in  the  preceding  sections. 

1.  The  simulations  performed  for  large  values  of  NA 
reproduce  the  work  performed  by  Pas  [Ref.  5]  and  therefore 
substantiate  those  results. 

2.  For  more  dilute  arrays,  because  the  problem  of 
particles  along  the  border  was  not  recognized,  Pas's 
calculations  for  first  NND  was  overestimated. 

3.  As  the  number  of  particles  utilized  in  any  simulation 
increases,  the  relative  number  of  particles  appearing  along 
the  border  decreases.  For  this  reason  the  errors  resulting  in 
calculations  of  first  NND  can  be  reduced  by  using  a  sufficient 
number  of  particles  (at  least  1000) .  Calculations  of  first 
NND  for  points  approached  the  theoretical  Poisson  distribution 
as  the  total  number  of  points  utilized  was  increased.  The 
same  can  be  said  for  the  distribution  of  finite  sized 
particles.  However  in  this  case,  as  the  total  number  of 
particles  increases  (approaches  1000  and  beyond) ,  6 
asymptotically  approaches  a  constant  value  which  lies  between 
the  minimum  possible  spacing  (2  x  rp)  and  the  spacing  for  a 
hexagonal  array.  Based  on  these  observations  it  is  possible 
to  compensate  for  the  statistical  inaccuracies  by  utilizing  at 
least  1000  particles  in  the  simulations. 


83 


4.  Attempts  to  compensate  for  errors  arising  in  the  first 
NND  calculations  when  using  a  small  number  of  particles  by 
eliminating  particles  along  the  border  was  only  partially 
successful.  Any  gains  made  using  the  procedure  to  reduce  the 
border  particle  problem  for  small  number  of  particles  was 
overtaken  by  simply  increasing  the  total  number  of  particles. 

5.  In  all  four  cases  analyzed  (random 
distribution/constant  size  particles,  non-random 
distribution/ constant  size  particles,  random 
distribution/ lognormal  particle  size  distribution,  and  non- 
random  distribution/ lognormal  particle  size  distribution)  the 
non-overlap  criterion  caused  the  distribution  to  deviate  from 
the  Poisson  to  that  of  a  hexagonal  array  as  the  density  of  the 
particles  increased.  Using  a  nominal  value  of  5000  particles 
the  model  very  closely  simulated  the  Poisson  distribution  in 
all  cases  for  dilute  arrays. 

6.  For  constant  size  particles  in  high  area  fractions  (af 
>  0.10),  for  both  random  and  non-random  distributions,  the 
most  frequently  occurring  first  NND  was  the  minimum  first  NND. 
This  was  due  to  the  non-overlap  criterion  and  that  the  minimum 
first  NND  can  not  be  less  than  the  particle  diameter. 

7.  The  non-random  distribution  of  points  and  constant 
size  particles  was  shifted  to  the  left  of  the  corresponding 
Poisson  distribution.  This  was  due  to  the  presence  of  the 
high  density  zones.   For  a  given  NA,  the  closer  packing  of  the 
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points  or  particles  in  the  high  density  zones  causes  NND  to 
shift  towards  lower  values.  This  was  also  reflected  in  a 
greater  frequency  of  the  minimum  first  NND  when  compared  to 
that  of  the  random  case. 

8.  The  random  distribution  of  particles  with  a  lognormal 
size  distribution  maintained  the  shape  of  a  Poisson  field  for 
all  area  fractions.  Unlike  the  case  for  constant  size 
particles,  the  minimum  first  NND  can  be  smaller  than  the  mean 
particle  diameter.  A  smaller  particle  can  always  fit  between 
two  or  more  larger  ones  thus  lowering  the  minimum  first  NND. 
This  also  resulted  in  a  more  constant  peak  frequency  of  first 
NND  for  all  ranges  of  area  fractions  examined.  These 
observations  held  true  as  well  for  the  case  of  non- 
random/  lognormal  particle  size  distribution. 

After  reviewing  all  the  accumulated  data  there  did  not 
appear  to  be  as  significant  a  difference  between  the  plotted 
distributions  of  the  random  and  non-random  cases  as  there  was 
between  the  plots  of  the  particles  themselves.  Based  then  on 
the  statistical  nature  of  this  work  it  appears  that  the  best 
dif ferentiater  between  random  and  non-random  particle 
distributions  is  the  human  eye. 

To  follow  up  on  this  work  the  following  studies  are 
recommended.  A  similar  analysis  should  be  made  for  the  three 
dimensional  case  in  order  to  parallel  the  actual  particle 
distributions.  Other  types  of  non-randomness,  such  as 
clusters,  should  be  explored  to  determine  their  effects  on 
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distributions.  Area  fractions  greater  than  0.50  should  be 
examined  as  well.  Finally,  an  analysis  should  be  made,  based 
on  results  of  this  work,  to  determine  the  feasibility  of 
developing  a  computer  program  which  could  determine  the  degree 
of  randomness  in  the  distributions  of  second  phase  particles. 
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APPENDIX  A.   PROGRAM  RANPART1 . FOR 


PROGRAM  RANPART1 

>f  variables: 
3.14159 
X  coordinate  of  an  interior  particle 

Y  coordinate  of  an  interior  particle 
X  coordinate  of  a  particle  along  the 
border 

Y  coordinate  of  a  particle  along  the 
border 

Area  percentage  of  the  particles 
Total  area  fraction  of  the  particles 
Min.  spacing  between  particles 
Particle  radius 

Defined  area  where  particles  are  positioned 
Area  occupied  by  each  particle 
Distance  between  interior  particles 
Distance  between  interior  particles 

Distance  between  an  interior  particle  and  a 
border  particle 
Nearest  neighbor  spacing 
Sum  of  nearest  neighbor  spacings 
Particle  diameter  plus  min  spacing 
Average  nearest  neighbor  spacing 
Horizontal  dimension  of  defined  area 
Vertical  dimension  of  defined  area 
Minimum  nearest  neighbor  spacing 
Distance  between  an  interior  particle  and 
its  nearest  border  particle 
Maximum  nearest  neighbor  spacing 
Size  of  each  range  between  dpmin  and  dpmax 
based  on  the  Sturgis  Rule 
Minimum  value  of  each  range  in  dpinc 
Number  of  particles  corresponding  to  each 
Delta 

Number  of  particles  per  unit  area 
Total  number  of  particles  generated 
Either  "A"  for  area  or  "p"  for  points, 
depending  on 
type  of  particles 

Either  "n"  for  a  preset  number  of  particles 
or  "r"  for  a  preset  particle  radius 

REAL  Pi,  Xpos(5200),  Ypos(5200),  XposB(250),  YposB(250)/ 
TPA,  AP,  Pspace,  Pradius,  TOTarea,  Parea,  dp,  dp2 ,  dp2B, 
dp2Bmin,  dp3(5200),  dp4 ,  dist,  dpavg,  X,  Y,  dpmin,  NPA, 
dpmax,  dpinc,  g,  Delta (30),  Np(30) 


* 

Glossary 

* 

Pi 

• 

Xpos 

* 

Ypos 

* 

XposB 

* 

* 

YposB 

* 

* 

AP 

* 

TPA 

* 

Pspace 

* 

Pradius 

4c 

TOTarea 

* 

Parea 

* 

dp 

* 

dp2 

* 

dp2B 

* 

* 

dp3 

* 

dp  4 

* 

dist 

* 

dpavg 

* 

X 

* 

Y 

* 

dpmin 

* 

dp2Bmin 

* 

* 

dpmax 

* 

dpinc 

* 

* 

Delta 

* 

Np 

* 

* 

NPA 

* 

Pnum 

* 

Ans 

* 

* 

* 

Part 

* 

* 
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* 

INTEGER  a,  b,  c,  d,  i,  j,  k,  1,  m,  n,  Pnum,  f,  h,  e 


CHARACTER  Ans*l,  Part*l 

PARAMETER  (Pi=3. 14159) 

DATA  n,  dp4,  e,  dpmin,  dpmax,  Np(l),  a,  b  / 
:1,  0.0,  0,  10.0,  0.001,  0,  1,  1  / 
* 

OPEN  (10, f ile='rplx.dat' , status=/new/ ) 

OPEN  (20,file='rply.daf , status='new' ) 

OPEN  (30,file='Dell.dat' , status='new' ) 

OPEN  (4  0,file='Npl.daf , status='new' ) 

OPEN  (50,file=/rpl.dat/ , status='new' ) 
* 

PRINT*,'  Enter  the  area  percentage  of  particles.' 

READ*,  AP 

PRINT*,'  Enter  the  length  of  the  X  -  scale.' 

READ*,  X 

PRINT*,'  Enter  the  length  of  the  Y  -  scale.' 

READ*,  Y 

TOTarea  =  X*Y 

TPA  =  AP/100*TOTarea 

PRINT*,'  Enter  the  minimum  spacing  between  particles.' 

READ*,  Pspace 

PRINT*,'  If  particle  generation  is  to  based  on  a  preset 
:  number  of 

PRINT*,'  particles  enter  "n" .   OR  If  particle 
: generation  is  to' 

PRINT*,'  be  based  on  a  preset  particle  radius  enter 

.     Ilv-ll  / 

•        ±.         m 

READ*, Part 

IF    (Part    .EQ.     'n')    THEN 

PRINT*, 'Enter  total  number  of  particles.' 
READ*, Pnum 
Parea  =  TPA/Pnum 
Pradius  =  SQRT(Parea/Pi) 
ELSEIF  (Part  .EQ.  'r')  THEN 

PRINT*,'  Enter  the  radius  of  the  particles.' 
READ*,  Pradius 
Parea  =  Pi* (Pradius**2) 
Pnum  =  NINT(TPA/Parea) 
ENDIF 

PRINT*,'  Are  the  particles  finite  areas  or  points? 
: Enter  "A"  or  "p" . ' 
READ* ,  Ans 
* 

*     Determine  random  position  of  particles. 


* 


PRINT*,'  Determining  random  position  of  particles.' 
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dist  =  2*Pradius  +  Pspace 
Xpos(l)  =  RRAND() 
DO  100  i  =  1,  Pnum 
150    CONTINUE 

Xpos(i)  =  RND()*X 
Ypos(i)  =  RND()*Y 
* 

*  Check  for  overlap  of  particles  with  a  finite  area 
* 

IF  (Ans  .EQ.  'p')  GOTO  250 
DO  200  k  =  1,  i-1 

dp  =  SQRT( (Xpos(i)  -  Xpos(k))**2  + 
:  (Ypos(i)  -  Ypos(k))**2) 

IF  (dp  .LT.  dist)  GOTO  150 
200        CONTINUE 
250      CONTINUE 
100    CONTINUE 

*  Segregate  particles  along  the  perimeter  from  those  in 

*  the  interior  of  the  specified  area. 
* 

PRINT*, '  Segregating  particles  on  the  border  from  the 
: interior. ' 
DO  27  5  d  =  1,  Pnum 

IF  ((Xpos(d)  .LE.  Pradius  +  Pspace)  .OR. 

(Xpos(d)  .GE.  X  -  Pradius  -  Pspace)  .OR. 
(Ypos(d)  .LE.  Pradius  +  Pspace)  .OR. 
(Ypos(d)  .GE.  Y  -  Pradius  -  Pspace))  THEN 
XposB(b)  =  Xpos(d) 
YposB(b)  =  Ypos(d) 
b  =  b  +  1 
ELSE 

Xpos(a)  =  Xpos(d) 
Ypos(a)  =  Ypos(d) 
WRITE (10,*)  Xpos(a) 
WRITE(20,*)  Ypos(a) 
a  =  a  +  1 
ENDIF 
275    CONTINUE 
* 

*  Determine  the  average  distance  between  nearest 

*  neighbors. 
* 

PRINT*,'  Determining  nearest  neighbor  distances.' 
DO  300  j  =  1,  Pnum  -  (b  -  1) 
dp3(n)  =  SQRT(X**2  +  Y**2) 
dp2Bmin  =  SQRT(X**2  +  Y**2) 
DO  350  c  =  1,  (b  -  1) 

dp2B  =  SQRT( (Xpos(j)  -  XposB(c))**2  + 
:  (Ypos(j)  -  YposB(c) )**2) 

dp2Bmin  =  MIN(dp2Bmin,  dp2B) 
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350 


400 


300 


500 

* 
* 
* 
* 

* 


600 


800 


700 


CONTINUE 

DO  400  1=1,  Pnum  -  (b  -  1) 

dp2  =  SQRT( (Xpos(l)  -  Xpos(j))**2  + 

(Ypos(l)  -  Ypos(j))**2) 
IF  (dp2  .NE.  0.0)  THEN 

dp3(n)  =  MIN(dp2,  dp3(n),  dp2Bmin) 
ENDIF 
CONTINUE 
n  =  n  +  1 
CONTINUE 
DO  500  m  =  1, 
dp4  =  dp4  + 
CONTINUE 
dpavg  =  dp4/(Pnum  -  (b  -  1)) 


Pnum  -  (b  -  1) 
dp3(m) 


Determine  the  number  of  particles  falling  in  each  range 
of  nearest  neighbor  spacings  (based  on  the  Sturgis 
rule) . 

PRINT*,'  Grouping  nearest  neighbor  values  into  ranges.' 
DO  600  f  =  1,  Pnum  -  1  -  (b  -  1) 

dpmin  =  MIN(dp3(f),  dp3(f+l),  dpmin) 

dpmax  =  MAX(dp3(f),  dp3(f+l),  dpmax) 
CONTINUE 

dpinc  =  (dpmax  -  dpmin) /(l  +  3 .  3*LOG10 (Pnum  -  (b  -  1))) 
DO  700  g  =  dpmin,  dpmax,  dpinc 

e  =  e  +  1 

Pnum  -  (b  -  1) 

.GE.  g)  .AND.  (dp3 (h)  .LT.  g  +  dpinc)) 


DO  800  h  =  1, 
IF  ((dP3(h) 
THEN 

Np(e) 
ENDIF 

CONTINUE 

Delta (e) 


=  Np(e)  +  1 


=  (2*g  +  dpinc) /2 
WRITE (30,*) Delta ( e ) 
WRITE(40,*)Np(e)/(Pnum  -  (b  -  1)) 
CONTINUE 

NPA  =  (Pnum  -  (b  -  1))/(X  -  2*Pradius  -  2*Pspace)/ 
(Y  -  2*Pradius  -  2*Pspace) 


WRITE(50,10) 
WRITE(50,20) 
WRITE(50,20) 
WRITE(50,20) 
WRITE(50,20) 
WRITE(50,20) 
WRITE(50,20) 
WRITE(50,20) 


'The  area  percentage  of  particles  is',  AP 

'X  dimension  is',X 

'Y  dimension  is',Y 

'Min  spacing  between  particles  is',Pspace 

'Particle  radius  is',  Pradius 

'Particle  area  is',  Parea 

'Total  number  of  particles  is',  Pnum 

'Number  of  particles  along  the  perimeter 


:is' ,b  -  1 

WRITE (50, 20) 'Number  of  particles  in  area  interior  is', a 
:-  1 
WRITE (50, 20) 'Number  particles  per  unit  area',  NPA 
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10 
20 


WRITE (50, 20) 'Average  spacing  between  nearest  neighbors 

is' ,  dpavg 

WRITE(50, 20) 'Min  nearest  neighbor  spacing 

WRITE (50, 20) 'Max  nearest  neighbor  spacing 


PRINT*, 
PRINT*, 
PRINT*, 
PRINT*, 
PRINT*, 
PRINT*, 
PRINT*, 
PRINT*, 
PRINT*, 
PRINT*, 
PRINT*, 
: dpavg 
PRINT*, 
PRINT*, 
FORMAT 
FORMAT 
END 


The  area  percentage  of  particles 


is' 
is' 
is' 


, dpmin 
, dpmax 
,  AP 


X  dimension  is',X 

Y  dimension  is',Y 

Min  spacing  between  particles  is',Pspace 

Particle  radius  is',  Pradius 

Particle  area  is',  Parea 

Total  number  of  particles  is',  Pnum 

Number  of  particles  along  the  perimeter', 

Number  of  particles  in  area  interior' ,  a 

Number  particles  per  unit  area 

Average  spacing  between  nearest 


is', 


NPA 
neighbors 


b  - 

-  1 

is1 


Min  nearest  neighbor  spacing 
Max  nearest  neighbor  spacing 
////// ,T8,A,F10. 4) 


is' , dpmin 
is' , dpmax 


(T8,A,F10.4) 
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APPENDIX  B.   PROGRAM  RANPART2.FOR 

PROGRAM  RANPART2 

of  variables: 
3.14159 

X  coordinate  of  an  interior  particle 
Y  coordinate  of  an  interior  particle 

X  coordinate  of  a  particle  along  the  border 

Y  coordinate  of  a  particle  along  the  border 
Area  percentage  of  the  particles 
Total  area  fraction  of  the  particles 
Min.  spacing  between  particles 
Particle  radius 
Defined  area  of  particles 
Area  occupied  by  each  particle 
Distance  between  interior  particles 
Distance  between  interior  particles 

Distance  between  an  interior  particle  and  a 
border  particle 
Nearest  neighbor  spacing 
Sum  of  nearest  neighbor  spacings 
Particle  diameter  plus  min  spacing 
Average  nearest  neighbor  spacing 
Horizontal  dimension  of  total  area 
Vertical  dimension  of  total  area 
Minimum  nearest  neighbor  spacing 
Distance  between  an  interior  particle  and 
its  nearest  border  particle 
Maximum  nearest  neighbor  spacing 
Size  of  each  range  between  dpmin  and  dpmax 
based  on  the  Sturgis  Rule 
Minimum  value  of  each  range  in  dpinc 
Number  of  particles  corresponding  to  each 
Delta 

Number  of  particles  per  unit  area 
Total  number  of  particles  generated 
Either  "A"  for  area  or  "p"  for  points, 
depending  on  type  of  particles 

Either  "n"  for  a  preset  number  of  particles 
or  "r"  for  a  preset  particle  radius 
Area  percentage  of  particles  in  each  of  the 
five  sub-areas 
Subtotal  of  particles  in  each  sub-area 

REAL  Pi,  Xpos(5200)/  Ypos(5200),  XposB(250),  YposB(250), 
TPA,  AP,  Pspace,  Pradius,  TOTarea,  Parea,  dp,  dp2 ,  dp2B, 
dp2Bmin,  dp3(5200),  dp4 ,  dist,  dpavg,  X,  Y,  dpmin,  NPA, 


* 

* 

Glossary 

* 

Pi 

* 

Xpos 

* 

Ypos 

* 

XposB 

* 

YposB 

* 

AP 

* 

TPA 

* 

Pspace 

* 

Pradius 

* 

TOTarea 

* 

Parea 

* 

dp 

* 

dp2 

* 

dp2B 

• 

* 

dp  3 

* 

dp4 

* 

dist 

* 

dpavg 

* 

X 

* 

Y 

* 

dpmin 

* 

dp2Bmin 

* 

* 

dpmax 

* 

dpinc 

* 

* 

Delta 

* 

Np 

* 

* 

NPA 

* 

Pnum 

* 

Ans 

* 

* 

Part 

* 

* 

AF1-5 

* 

* 

Al-5 

* 
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* 

* 

* 
* 


:dpmax,  dpinc,  g,  Delta(30),  Np(30),  Al,  AF1,  A2 ,  AF2 , 
:A3,  AF3,  A4 ,  AF4 ,  A5,  AF5 

INTEGER  a,  b,  c,  d,  i,  j,  k,  1,  m,  n,  Pnum,  f,  h,  e 

CHARACTER  Ans*l,  Part*l 

PARAMETER  (Pi=3. 14159) 

DATA  n,  dp4 ,  e,  dpmin,  dpmax,  Np(l),  a,  b,  Al,  A2 ,  A3, 
:A4,  A5  / 

:1,  0.0,  0,  10.0,  0.001,  0,  1,  1,  0.0,  0.0,  0.0,  0.0, 
:0.0/ 

OPEN  (10,f ile=' rp2x.dat' , status='new' ) 

OPEN  (20,file='rp2y.daf , status='new' ) 

OPEN  (30,file='Del2.daf , status='new' ) 

OPEN  (40,file='Np2.dat' , status='new' ) 

OPEN  (50,file='rp2.daf , status='new' ) 

PRINT*,'  Enter  the  area  percentage  of  particles.' 

READ*,  AP 

PRINT*,'  Enter  5  sub-area  percentages  of  particles.' 

READ*,  AF1,  AF2,  AF3 ,  AF4 ,  AF5 

AF1  =  AF1/100 

AF2  =  AF2/100 

AF3  =  AF3/100 

AF4  =  AF4/100 

AF5  =  AF5/100 

PRINT*,'  Enter  the  length  of  the  X  -  scale.' 

READ*,  X 

PRINT*,'  Enter  the  length  of  the  Y  -  scale.' 

READ*,  Y 

TOTarea  =  X*Y 

TPA  =  AP/100*TOTarea 

PRINT*,'  Enter  the  minimum  spacing  between  particles.' 

READ* ,  Pspace 

PRINT*,'  If  particle  generation  is  to  based  on  a  preset 
: number  of ' 

PRINT*,'  particles  enter  "n" .   OR  If  particle 
: generation  is  to' 

PRINT*,'  be  based  on  a  preset  particle  radius  enter 

.  My  II    / 

READ*, Part 

IF  (Part  .EQ.  'n')  THEN 

PRINT*, 'Enter  total  number  of  particles.' 

READ* , Pnum 

Parea  =  TPA/ Pnum 

Pradius  =  SQRT(Parea/Pi) 
ELSEIF  (Part  .EQ.  'r')  THEN 

PRINT*,'  Enter  the  radius  of  the  particles.' 
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READ*,  Pradius 

Parea  =  Pi* (Pradius**2) 

Pnum  =  NINT(TPA/Parea) 
ENDIF 

PRINT*,'  Are  the  particles  finite  areas  or  points? 
: Enter  "A"  or  "p". ' 
READ* ,  Ans 

*  Determine  random  position  of  particles. 
* 

PRINT*,'  Determining  random-banded  position  of 

: particles. ' 

dist  =  2*Pradius  +  Pspace 

Xpos(l)  =  RRAND() 

DO  100  i  =  1,  Pnum 

150        CONTINUE 

Xpos(i)  =  RND()*X 

Ypos(i)  =  RND()*Y 
* 

*  Check  for  overlap  of  particles  with  a  finite  area 
* 

IF  (Ans  .EQ.  'p')  GOTO  250 
DO  200  k  =  1,  i  -  1 

dp  =  SQRT( (Xpos(i)  -  Xpos(k))**2  + 

(Ypos(i)  -  Ypos(k))**2) 
IF  (dp  .LT.  dist)  GOTO  150 
2  00        CONTINUE 
2  50    CONTINUE 

IF  ((Ypos(i)  .GE.  0.0)  .AND.  (Ypos(i)  .LT.  0.2*Y))  THEN 
IF  (Al  .GT.  TPA*AF1)  GOTO  150 
Al  =  Al  +  Parea 
ELSEIF  ((Ypos(i)  .GE.  0.2*Y)  .AND.  (Ypos(i)  .LT.  0.4*Y)) 
:THEN 

IF  (A2  .GT.  TPA*AF2)  GOTO  150 
A2  =  A2  +  Parea 
ELSEIF  ((Ypos(i)  .GE.  0.4*Y)  .AND.  (Ypos(i)  .LT.  0.6*Y)) 
:THEN 

IF  (A3  .GT.  TPA*AF3)  GOTO  150 
A3  =  A3  +  Parea 
ELSEIF  ((Ypos(i)  .GE.  0.6*Y)  .AND.  (Ypos(i)  .LT.  0.8*Y)) 
:THEN 

IF  (A4  .GT.  TPA*AF4)  GOTO  150 
A4  =  A4  +  Parea 
ELSEIF  ((Ypos(i)  .GE.  0.8)  .AND.  (Ypos(i)  .LE.  Y) )  THEN 
IF  (A5  .GT.  TPA*AF5)  GOTO  150 
A5  =  A5  +  Parea 
ENDIF 
100    CONTINUE 
* 

*  Segregate  particles  along  the  perimeter  from  those  in 

*  the  interior  of  the  specified  area. 
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PRINT*, '  Segregating  particles  on  the  border  from  the 

interior. ' 

DO  2  75  d  =  1,  Pnum 

IF  ((Xpos(d)  .LE.  Pradius  +  Pspace)  .OR. 

(Xpos(d)  .GE.  X  -  Pradius  -  Pspace)  .OR. 
(Ypos(d)  .LE.  Pradius  +  Pspace)  .OR. 
(Ypos(d)  .GE.  Y  -  Pradius  -  Pspace))  THEN 
XposB(b)  =  Xpos(d) 
YposB(b)  =  Ypos(d) 
b  =  b  +  1 
ELSE 

Xpos(a)  =  Xpos(d) 
Ypos(a)  =  Ypos(d) 
WRITE(10,*)  Xpos(a) 
WRITE(20,*)  Ypos(a) 
a  =  a  +  1 
END  IF 
CONTINUE 

Determine  the  average  distance  between  nearest 
neighbors. 

PRINT*,'  Determining  nearest  neighbor  distances.' 
DO  300  j  =  1,  Pnum  -  (b  -  1) 
dp3(n)  =  SQRT(X**2  +  Y**2) 
dp2Bmin  =  SQRT(X**2  +  Y**2) 
DO  350  C  =  1,  (b  -  1) 

dp2B  =  SQRT( (Xpos(j)  -  XposB(c))**2  + 

(Ypos(j)  -  YposB(c) )**2) 
dp2Bmin  =  MIN(dp2Bmin,  dp2B) 
CONTINUE 
DO  400  1=1,  Pnum  -  (b  -  1) 

dp2  =  SQRT( (Xpos(l)  -  Xpos(j))**2  + 

(Ypos(l)  -  Ypos(j) )**2) 
IF  (dp2  .NE.  0.0)  THEN 

dp3(n)  =  MIN(dp2,  dp3 (n) ,  dp2Bmin) 
ENDIF 
CONTINUE 
n  =  n  +  1 
CONTINUE 

DO  500  m  =  1,  Pnum  -  (b  -  1) 
dp4  =  dp4  +  dp3  (m) 
)0    CONTINUE 

dpavg  =  dp4/(Pnum  -  (b  -  1)) 
* 

*  Determine  the  number  of  particles  falling  in  each  range 

*  of  nearest  neighbor  spacings  (based  on  the  Sturgis 

*  rule) . 

k 

DO  600  f  =  1,  Pnum  -  1  -  (b  -  1) 
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600 


800 


700 


dpmin  =  MIN(dp3(f),  dp3(f+l),  dpmin) 
dpmax  =  MAX(dp3(f),  dp3(f+l),  dpmax) 

CONTINUE 

dpinc  =  (dpmax  -  dpmin) /(l  +  3 . 3*LOG10 (Pnum  - 

DO  700  g  =  dpmin,  dpmax,  dpinc 
e  =  e  +  1 


(b  -  1))) 


DO  800  h 

=  1, 

Pnum 

- 

(b  -  1) 

IF  ((dp3(h) 

.GE. 

g) 

.AND. 

THEN 

Np(e) 

=  Np( 

e)  + 

i 

ENDIF 

CONTINUE 

Delta (e) 

=  (2* 

g  +  < 

ipinc) /2 

(dp3(h)  .LT.  g  +  dpinc)) 


(b  -  1)) 


is 


WRITE(30,*)Delta(e) 

WRITE(40,*)Np(e)/(Pnum  - 
CONTINUE 
NPA  =  (Pnum  -  (b  -  1))/(X  -  2*Pradius  -  2*Pspace)/ 

(Y  -  2*Pradius  -  2*Pspace) 
WRITE(50,10 
WRITE(50,20 
WRITE(50,20 
WRITE(50,20 
WRITE(50,20 
WRITE(50,20 
WRITE(50,20 
WRITE(50,20 
is' ,b  -  1 
WRITE(50,20 
:-  1 

WRITE(50,20 
WRITE(50,20 
is' ,  dpavg 
WRITE(50,20 
WRITE(50,20 
PRINT*, 
PRINT*, 
PRINT* , 
PRINT*, 
PRINT*, 
PRINT*, 
PRINT*, 
PRINT* , 
1 

PRINT*, 
PRINT*, 
PRINT*, 
dpavg 
PRINT*, 
PRINT*, 


AP 


The  area  percentage  of  particles 

X  dimension  is',X 

Y  dimension  is',Y 

Min  spacing  between  particles  is',Pspace 

Particle  radius  is',  Pradius 

Particle  area  is',  Parea 

Total  number  of  particles  is' ,  Pnum 

Number  of  particles  along  the  perimeter 

Number  of  particles  in  area  interior  is', a 

Number  particles  per  unit  area',  NPA 
Average  spacing  between  nearest  neighbors 

Min  nearest  neighbor  spacing  is', dpmin 
Max  nearest  neighbor  spacing  is', dpmax 


is 


AP 


The  area  percentage  of  particles 

X  dimension  is',X 

Y  dimension  is',Y 

Min  spacing  between  particles  is',Pspace 


Particle 
Particle 


radius  is',  Pradius 


area  is 


Parea 


Pnum 


Total  number  of  particles  is 

Number  of  particles  along  the  perimeter',  b  - 


Number 
Number 


of  particles  in  area  interior',  a  -  1 
particles  per  unit  area  is' 


NPA 


Average  spacing  between  nearest  neighbors  is', 


Min  nearest  neighbor  spacing  is', dpmin 
Max  nearest  neighbor  spacing  is', dpmax 
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10     FORMAT  (  / / / / / / ,T8,A,F10.4) 

20     FORMAT  (T8,A,F10.4) 
END 
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APPENDIX  C.   PROGRAM  RANPART3.FOR 


PROGRAM  RANPART3 


of  variables: 
3.14159 

X  coordinate  of  an  interior  particle 
Y  coordinate  of  an  interior  particle 
X  coordinate  of  a  particle  along  the  border 
Y  coordinate  of  a  particle  along  the  border 
Area  percentage  of  the  particles 
Total  area  occupied  by  the  particles 
Min.  spacing  between  particles 
Defined  area  of  particles 
Area  occupied  by  each  particle 
Distance  between  interior  particles 
Distance  between  interior  particles 
Distance  between  an  interior  particle  and  a 
border  particle 
Nearest  neighbor  spacing 
Sum  of  nearest  neighbor  spacings 
Particle  diameter  plus  min  spacing 
Average  nearest  neighbor  spacing 
Horizontal  dimension  of  total  area 
Vertical  dimension  of  total  area 
Minimum  nearest  neighbor  spacing 
Distance  between  an  interior  particle  and 
its  nearest  border  particle 
Maximum  nearest  neighbor  spacing 
Size  of  each  range  between  dpmin  and  dpmax 
based  on  the  Sturgis  Rule 
Minimum  value  of  each  range  in  dpinc 
Number  of  particles  corrosponding  to  each 
Delta 

Number  of  particles  per  unit  area 
Total  number  of  particles  generated 
Generated  particle  radius  from  lonormal 
distribution 
Stored  particle  radius 
Stored  border  particle  radius 
Standard  deviation 
Log  mean  particle  radius 
Maximum  particle  radius 
Minimum  particle  radius 
Largest  border  particle 

REAL  Pi,  Xpos(lOOOO),  Ypos(lOOOO),  Rad(lOOOO), 
XposB(300),  YposB(300),  RadB(300),  R,  TPA,  AP,  Pspace, 


* 

Glossary 

* 

Pi 

* 

Xpos 

* 

Ypos 

* 

XposB 

* 

YposB 

* 

AP 

* 

TPA 

* 

Pspace 

* 

TOTarea 

* 

Parea 

* 

dp 

* 

dp2 

* 

dp2B 

* 

* 

dp  3 

* 

dp4 

* 

dist 

* 

dpavg 

* 

X 

* 

Y 

* 

dpmin 

* 

dp2Bmin 

* 

* 

dpmax 

* 

dpinc 

* 

* 

Delta 

* 

Np 

* 

* 

NPA 

* 

Pnum 

* 

R 

* 

* 

Rad ( i ) 

* 

RadB(i) 

* 

S 

* 

u 

* 

Radmax 

* 

Radmin 

* 

RadBmax 

* 
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:TOTarea,  Parea,  dp,  TRad,  dp2 ,  dp2B,  dp2Bmin,  u,  S, 
:dp3(10000),  dp4,  dist,  dpavg,  X,  Y,  dpmin,  NPA, 
:dpmax,  dpinc,  g,  Delta(50),  Np(50),  Pnum,  RadBmax, 
:Radmax,  Radmin 

INTEGER  a,  b,  c,  d,  e,  f,  h,  i,  j,  k,  1,  m,  n,  p,  q,  nq 

PARAMETER  (Pi=3. 14159) 

DATA  n,  dp4,  e,  dpmin,  dpmax,  Np(l),  a,  b,  Parea,  TRad, 
:RadBmax,  Radmax,  Radmin,  nr,  nq  / 

:1,  0.0,  0,  10.0,  0.001,  0,  1,  1,  0.0,  0.0,  .00001, 
: .00001,  1000,  0,0/ 

OPEN  (10,file='rp3x.dat' , status='NEW ) 

OPEN  (20, file=' rp3y.dat' , status='NEW ) 

OPEN  (30,file=/Del3.dat' , status=/NEW/ ) 

OPEN  (40,file=,Np3.dat/ , status=,NEW/ ) 

OPEN  (50, file=' rp3.dat' , status=,NEW ) 

OPEN  (60^116=' Rad3.dat'  ,  status='NEW  ) 

PRINT*,'  Enter  the  area  percentage  of  particles.' 
READ*,  AP 

PRINT*,'  Enter  the  length  of  the  X  -  scale.' 
READ* ,  X 

PRINT*,'  Enter  the  length  of  the  Y  -  scale.' 
READ*,  Y 
TOTarea  =  X*Y 
TPA  =  AP/100*TOTarea 

PRINT*,'  Enter  the  mean  particle  radius  and  standard 
: deviation. . . ' 
READ*,  u,  S 

PRINT*,'  Enter  the  minimum  spacing  between  particles.' 
READ*,  Pspace 

PRINT*,'  Enter  the  total  number  of  particles...' 
READ* , Pnum 

Determine  random  position  and  radius  of  each  particle. 

ISEED  =  0 

CALL  RNSET  (ISEED) 

DO  100  i  =  1,  Pnum 

CONTINUE 

q  =  0 

CALL  RNLNL  (1,  u,  S,  R) 

CONTINUE 

q  =  q  +  1 

IF  (q  .GE.  600)  THEN 

nq  =  nq  +  1 

GOTO  175 

END  IF 
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Xpos(i)  =  RNUNF()*X 
Ypos(i)  =  RNUNF()*Y 
Rad(i)  =  R 
* 

*  Check  for  overlap  of  particles 
* 

DO  200  k  =  1,  i  -  1 

dp  =  SQRT( (Xpos(i)  -  Xpos(k))**2  + 
:  (Ypos(i)  -  Ypos(k))**2) 

dist  =  Rad(i)  +  Rad(k)  +  Pspace 
IF  (dp  .LT.  dist)  GOTO  150 
2  00         CONTINUE 

Parea  =  Parea  +  Pi*Rad(i)**2 
IF  (Parea  .GE.  TPA)  THEN 
Pnum  =  i 
GOTO  2  50 
ENDIF 
100    CONTINUE 
2  50    CONTINUE 
* 

*  Segregate  particles  along  the  perimeter  from  those  in 

*  the  interior  of  the  specified  area. 
* 

DO  275  d  =  1,  Pnum 

IF  ((Xpos(d)  .LE.  Rad(d)  +  Pspace)  .OR. 

(Xpos(d)  .GE.  X  -  Rad(d)  -  Pspace)  .OR. 
(Ypos(d)  .LE.  Rad(d)  +  Pspace)  .OR. 
(Ypos(d)  .GE.  Y  -  Rad(d)  -  Pspace))  THEN 
RadBmax  =  MAX(Rad(d),  RadBmax) 
ENDIF 
275    CONTINUE 

DO  280  p  =  1,  Pnum 

IF  ((Xpos(p)  .LE.  RadBmax  +  Pspace)  .OR. 

(Xpos(p)  .GE.  X  -  RadBmax  -  Pspace)  .OR. 
(Ypos(p)  .LE.  RadBmax  +  Pspace)  .OR. 
(Ypos(p)  .GE.  Y  -  RadBmax  -  Pspace))  THEN 
XposB(b)  =  Xpos(p) 
YposB(b)  =  Ypos(p) 
RadB(b)  =  Rad(p) 
b  =  b  +  1 
ELSE 

Xpos(a)  =  Xpos(p) 
Ypos(a)  =  Ypos(p) 
Rad(a)  =  Rad(p) 
TRad  =  TRad  +  Rad(a) 
WRITE(10,*)  Xpos(a) 
WRITE(20,*)  Ypos(a) 
WRITE(60,*)  Rad(a) 
a  =  a  +  1 
ENDIF 
2  80    CONTINUE 
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Determine  the  average  distance  between  nearest 
neighbors. 

DO  300  j  =  1,  Pnum  -  (b  -  1) 
dp3(n)  =  SQRT(X**2  +  Y**2) 
dp2Bmin  =  SQRT(X**2  +  Y**2) 
DO  350  c  =  1,  (b  -  1) 

dp2B  =  SQRT( (Xpos(j)  -  XposB(c))**2  + 
:  (Ypos(j)  -  YposB(c) )**2) 

dp2Bmin  =  MIN(dp2Bmin,  dp2B) 
CONTINUE 
DO  400  1=1,  Pnum  -  (b  -  1) 

dp2  =  SQRT( (Xpos(l)  -  Xpos(j))**2  + 
:  (Ypos(l)  -  Ypos(j) )**2) 

IF  (dp2  .NE.  0.0)  THEN 

dp3(n)  =  MIN(dp2,  dp3 (n) ,  dp2Bmin) 
ENDIF 
CONTINUE 
n  =  n  +  1 
CONTINUE 
DO  500  m  =  1,  Pnum  -  (b  -  1) 

dp4  =  dp4  +  dp3 (m) 
CONTINUE 
dpavg  =  dp4/(Pnum  -  (b  -  1)) 

Determine  the  number  of  particles  falling  in  each 
range  of  nearest  neighbor  spacings  (based  on  the  Sturgis 
rule) . 

DO  600  f  =  1,  Pnum  -  (b  -  1) 
dpmin  =  MIN(dp3(f),  dpmin) 
dpmax  =  MAX(dp3(f),  dpmax) 
Radmin  =  MIN(Rad(f) ,  Radmin) 
Radmax  =  MAX(Rad(f),  Radmax) 
CONTINUE 

dpinc  =  (dpmax  -  dpmin) /(l  +  3 . 3*LOG10 (Pnum  -  (b  -  1) ) ) 
DO  700  g  =  dpmin,  dpmax,  dpinc 
e  =  e  +  1 
DO  800  h  =  1,  Pnum  -  (b  -  1) 

IF  ((dp3(h)  .GE.  g)  .AND.  (dp3 (h)  .LT.  g  +  dpinc)) 
:     THEN 

Np(e)  =  Np(e)  +  1 
ENDIF 
CONTINUE 

Delta(e)  =  (2*g  +  dpinc) /2 
WRITE (30,*) Delta (e) 
WRITE(40,*)Np(e) /(Pnum  -  (b  -  1) ) 
CONTINUE 

NPA  =  (Pnum  -  (b  -  1) ) / (X  -  2*RadBmax  -  2*Pspace)/ 
:       (Y  -  2*RadBmax  -  2*Pspace) 
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10 
20 
30 


WRITE ( 
WRITE ( 
WRITE ( 
WRITE ( 
WRITE ( 

"  1)) 
WRITE ( 
WRITE ( 
WRITE ( 

-  1 
WRITE ( 

-  1 

WRITE ( 
WRITE ( 
WRITE ( 
WRITE ( 
WRITE ( 
WRITE ( 
PRINT* 
PRINT* 
PRINT* 
PRINT* 
PRINT* 
PRINT* 
PRINT* 
PRINT* 
PRINT* 
PRINT* 
PRINT* 
PRINT* 
PRINT* 
PRINT* 
PRINT* 
FORMAT 
FORMAT 
FORMAT 
END 


50,10 
50,20 
50,20 
50,20 
50,20 

50,20 
50,20 
50,30 

50,30 


The  area  percentage  of  particles  is',  AP 

X  dimension  is',X 

Y  dimension  is',Y 

Min  spacing  between  particles  is',Pspace 

Mean  particle  radius  is',  TRad/ (Pnum  -  (b 

Total  area  occupied  by  particles',  Parea 
Total  number  of  particles  is',  Pnum 
Number  of  particles  along  the  border  is',b 

Number  of  particles  in  area  interior  is', a 


50, 20) 'Number  particles  per  unit  area',  NPA 
50,20)  'Average  nearest  neighbor  spacing  is' ,  dpavg 
50,20) 'Min  nearest  neighbor  spacing  is',dpmin 
50,20) 'Max  nearest  neighbor  spacing  is',dpmax 
50,20) 'Max  particle  radius  is',Radmax 
50,20) 'Min  particle  radius  is',Radmin 

The  area  percentage  of  particles  is',  AP 
X  dimension  is',X 
Y  dimension  is',Y 

Min  spacing  between  particles  is',Pspace 
Mean  particle  radius  is',  TRad/ (Pnum  -  (b  -  1)  ) 
Total  area  occupied  by  particles  is',  Parea 
Total  number  of  particles  is',  Pnum 
Number  of  particles  along  the  border',  b  -  1 
Number  of  particles  in  area  interior' ,  a  -  1 
Number  particles  per  unit  area  is',  NPA 
Average  nearest  neighbor  spacing  is',  dpavg 
Min  nearest  neighbor  spacing  is',dpmin 
Max  nearest  neighbor  spacing  is',dpmax 
Max  particle  radius  is',Radmax 
Min  particle  radius  is',Radmin 
(T6,A,T50,F10.4) 
(T6,A,T50,F10.4) 
(T6,A,T46,I9) 
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APPENDIX  D.   PROGRAM  RANPART4.FOR 


PROGRAM  RANPART4 


* 

Glossary 

* 

Pi 

* 

Xpos 

* 

Ypos 

* 

XposB 

* 

YposB 

* 

AP 

* 

TPA 

* 

Pspace 

* 

TOTarea 

* 

Parea 

* 

dp 

* 

dp2 

* 

dp2B 

* 

* 

dp  3 

* 

dp4 

* 

dist 

* 

dpavg 

* 

X 

* 

Y 

* 

dpmin 

* 

dp2Bmin 

* 

* 

dpmax 

* 

dpinc 

* 

* 

Delta 

* 

Np 

* 

* 

NPA 

* 

Pnum 

* 

R 

* 

* 

Rad(i) 

* 

RadB(i) 

* 

S 

* 

u 

* 

Radmax 

* 

Radmin 

* 

RadBmax 

* 

AF1-5 

* 

* 

Al-5 

of  variables: 
3.14159 

X  coordinate  of  an  interior  particle 
Y  coordinate  of  an  interior  particle 
X  coordinate  of  a  particle  along  the  border 
Y  coordinate  of  a  particle  along  the  border 
Area  percentage  of  the  particles 
Total  area  occupied  by  the  particles 
Min.  spacing  between  particles 
Defined  area  of  particles 
Area  occupied  by  each  particle 
Distance  between  interior  particles 
Distance  between  interior  particles 
Distance  between  an  interior  particle  and  a 
border  particle 
Nearest  neighbor  spacing 
Sum  of  nearest  neighbor  spacings 
Particle  diameter  plus  min  spacing 
Average  nearest  neighbor  spacing 
Horizontal  dimension  of  total  area 
Vertical  dimension  of  total  area 
Minimum  nearest  neighbor  spacing 
Distance  between  an  interior  particle  and 
its  nearest  border  particle 
Maximum  nearest  neighbor  spacing 
Size  of  each  range  between  dpmin  and  dpmax 
based  on  the  Sturgis  Rule 
Minimum  value  of  each  range  in  dpinc 
Number  of  particles  corresponding  to  each 
Delta 

Number  of  particles  per  unit  area 
Total  number  of  particles  generated 
Generated  particle  radius  from  lonormal 
distribution 
Stored  particle  radius 
Stored  border  particle  radius 
Standard  deviation 
Log  mean  particle  radius 
Maximum  particle  radius 
Minimum  particle  radius 
Largest  border  particle 

Area  percentage  of  particles  in  each  of  the 
five  sub-areas 
Subtotal  of  particles  in  each  sub-area 
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** 


* 


REAL  Pi,  Xpos(lOOOO),  Ypos(lOOOO),  Rad(lOOOO), 
:XposB(300),  YposB(300),  RadB(300),  R,  TPA,  AP,  Pspace, 
:u,  TOTarea,  Parea,  dp,  TRad,  dp2 ,  dp2B,  dp2Bmin,  S, 
:dp3(10000),  dp4,  dist,  dpavg,  X,  Y,  dpmin,  NPA,  dpmax, 
:dpinc,  g,  Delta(50) ,  Np(50) ,  Pnum,  RadBmax,  Radmax, 
:Radmin,  Al,  AF1,  A2 ,  AF2 ,  A3,  AF3 ,  A4 ,  AF4 ,  A5 ,  AF5 

INTEGER  a,  b,  c,  d,  e,  f,  h,  i,  j,  k,  1,  m,  n,  p,  q,  nq 

PARAMETER    (Pi=3. 14159) 

DATA  n,    dp4,    e,    dpmin,    dpmax,    Np(l),    a,    b,    Parea,    TRad, 
: RadBmax,    Radmax,    Radmin,    Al,    A2,    A3,    A4 ,    A5,    nr,    nq   / 
:1,    0.0,    0,     10.0,    0.001,    0,     1,     1,    0.0,    0.0,     .00001, 
:. 00001,     1000,    0.0,    0.0,    0.0,    0.0,     0.0,    0,    0    / 

OPEN  (10,file='rp4x.daf , status=/NEW/ ) 

OPEN  (20,file='rp4y.daf , status='NEW ) 

OPEN  (30,file='Del4.daf , status='NEW ) 

OPEN  (40,file='Np4.daf , status='NEW ) 

OPEN  (50,file='rp4.dat',status='NEW) 

OPEN  (60,file='Rad4.daf , status='NEW ) 

PRINT*,'  Enter  the  area  percentage  of  particles.' 
READ*,  AP 

PRINT*,'  Enter  5  sub-area  percentages  of  particles' 
READ*,  AF1,  AF2,  AF3 ,  AF4 ,  AF5 
AF1  =  AF1/100 
AF2  =  AF2/100 
AF3  =  AF3/100 
AF4  =  AF4/100 
AF5  =  AF5/100 

PRINT*,'  Enter  the  length  of  the  X  -  scale.' 
READ* ,  X 

PRINT*,'  Enter  the  length  of  the  Y  -  scale.' 
READ*,  Y 
TOTarea  =  X*Y 
TPA  =  AP/100*TOTarea 

PRINT*,'  Enter  the  mean  particle  radius  and  standard 
: deviation. . . ' 
READ*,  u,  S 

PRINT*,'  Enter  the  minimum  spacing  between  particles.' 
READ* ,  Pspace 

PRINT*,'  Enter  the  total  number  of  particles...' 
READ* , Pnum 


*     Determine  random  position  and  radius  of  each  particle. 
* 

ISEED  =  0 
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CALL  RNSET  (ISEED) 
DO  100  i  =  1,  Pnum 
CONTINUE 

q  =  0 

CALL  RNLNL  (1,  u,  S,  R) 

CONTINUE 

q  =  q  +  1 

IF  (q  .GE.  600)  THEN 

nq  =  nq  +  1 

GOTO  175 

END  IF 

Xpos(i)  =  RNUNF()*X 

Ypos(i)  =  RNUNF()*Y 

Rad(i)  =  LOG(R) 

Check  for  overlap  of  particles 

DO  200  k  =  1,  i  -  1 

dp  =  SQRT( (Xpos(i)  -  Xpos(k))**2  + 
:  (Ypos(i)  -  Ypos(k))**2) 

dist  =  Rad(i)  +  Rad(k)  +  Pspace 
IF  (dp  .LT.  dist)  GOTO  150 
CONTINUE 
IF  ((Ypos(i)  .GE.  0.0)  .AND.  (Ypos(i)  .LT.  0.2*Y))  THEN 
IF  (Al  .GT.  TPA*AF1)  GOTO  150 
Al  =  Al  +  Pi*Rad(i)**2 
ELSEIF  ((Ypos(i)  .GE.  0.2*Y)  .AND.  (Ypos(i)  .LT.  0.4*Y)) 
:THEN 

IF  (A2  .GT.  TPA*AF2)  GOTO  150 
A2  =  A2  +  Pi*Rad(i)**2 
ELSEIF  ((Ypos(i)  .GE.  0.4*Y)  .AND.  (Ypos(i)  .LT.  0.6*Y)) 
:THEN 

IF  (A3  .GT.  TPA*AF3)  GOTO  150 
A3  =  A3  +  Pi*Rad(i)**2 

ELSEIF  ((Ypos(i)  .GE.  0.6*Y)  .AND.  (Ypos(i)  .LT.  0.8*Y)) 
:THEN 

IF  (A4  .GT.  TPA*AF4)  GOTO  150 
A4  =  A4  +  Pi*Rad(i)**2 
ELSEIF  ((Ypos(i)  .GE.  0.8*Y)  .AND.  (Ypos(i)  .LE.  Y) ) 
:THEN 

IF  (A5  .GT.  TPA*AF5)  GOTO  150 
A5  =  A5  +  Pi*Rad(i)**2 
END  IF 

Parea  =  Parea  +  Pi*Rad(i)**2 
IF  (Parea  .GE.  TPA)  THEN 
Pnum  =  i 
GOTO  250 
END  IF 


100    CONTINUE 
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250 

* 

* 

* 

* 


275 


280 

* 

* 


CONTINUE 

Segregate  particles  along  the  border  from  those  in 
the  in  of  the  specified  area. 


OR. 


THEN 


RadBmax  +  Pspace)  .OR. 
X  -  RadBmax  -  Pspace)  .OR. 
RadBmax  +  Pspace)  .OR. 
Y  -  RadBmax  -  Pspace) )  THEN 


350 


400 


DO  275  d  =  1,  Pnum 

IF  ((Xpos(d)  .LE.  Rad(d)  +  Pspace)  .OR. 
(Xpos(d)  .GE.  X  -  Rad(d)  -  Pspace) 
(Ypos(d)  .LE.  Rad(d)  +  Pspace)  .OR. 
(Ypos(d)  .GE.  Y  -  Rad(d)  -  Pspace)) 
RadBmax  =  MAX(Rad(d) ,  RadBmax) 
ENDIF 
CONTINUE 

DO  280  p  =  1,  Pnum 
IF  ((Xpos(p)  .LE 
(Xpos(p)  .GE. 
(Ypos(p)  .LE. 
(Ypos(p)  .GE. 
XposB(b)  =  Xpos(p) 
YposB(b)  =  Ypos(p) 
RadB(b)  =  Rad(p) 
b  =  b  +  1 
ELSE 

Xpos(a)  =  Xpos(p) 
Ypos(a)  =  Ypos(p) 
Rad(a)  =  Rad(p) 
TRad  =  TRad  +  Rad(a) 
WRITE(10,*)  Xpos(a) 
WRITE(20,*)  Ypos(a) 
WRITE (60,*)  Rad(a) 
a  =  a  +  1 
ENDIF 
CONTINUE 


Determine  the  average  distance  between  nearest 

DO  300  j  =  1,  Pnum  -  (b  -  1) 
dp3(n)  =  SQRT(X**2  +  Y**2) 
dp2Bmin  =  SQRT(X**2  +  Y**2) 
DO  350  c  =  1,  (b  -  1) 

dp2B  =  SQRT( (Xpos(j)  -  XposB(c))**2  + 

(Ypos(j)  -  YposB(c) )**2) 
dp2Bmin  =  MIN(dp2Bmin,  dp2B) 
CONTINUE 
DO  400  1=1,  Pnum  -  (b  -  1) 

dp2  =  SQRT( (Xpos(l)  -  Xpos(j))**2  + 

(Ypos(l)  -  Ypos(j) )**2) 
IF  (dp2  .NE.  0.0)  THEN 

dp3(n)  =  MIN(dp2,  dp3(n),  dp2Bmin) 
ENDIF 
CONTINUE 
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300 


500 

* 
* 


600 


800 


700 


n  =  n  +  1 
CONTINUE 
DO  500  m  =  1,  Pnum  -  (b  -  1) 

dp4  =  dp4  +  dp3 (m) 
CONTINUE 
dpavg  =  dp4/(Pnum  -  (b  -  1)) 

Determine  the  number  of  particles  falling  in  each  range 
of  nearest  neighbor  spacings  (based  on  the  Sturgis 
rule) . 

DO  600  f  =  1,  Pnum  -  (b  -  1) 
dpmin  =  MIN(dp3(f),  dpmin) 
dpmax  =  MAX(dp3(f),  dpmax) 
Radmin  =  MIN(Rad(f),  Radmin) 
Radmax  =  MAX(Rad(f),  Radmax) 
CONTINUE 

dpinc  =  (dpmax  -  dpmin) /(l  +  3 . 3*LOG10 (Pnum  -  (b  -  1) ) ) 
DO  7  00  g  =  dpmin,  dpmax,  dpinc 
e  =  e  +  1 
DO  800  h  =  1,  Pnum  -  (b  -  1) 

IF  ((dp3(h)  .GE.  g)  .AND.  (dp3 (h)  .LT.  g  +  dpinc)) 
:     THEN 

Np(e)  =  Np(e)  +  1 
END  IF 
CONTINUE 

Delta(e)  =  (2*g  +  dpinc) /2 
WRITE (30,*) Delta (e) 
WRITE(40,*)Np(e)/(Pnum  -  (b  -  1)  ) 
CONTINUE 

NPA  =  (Pnum  -  (b  -  1))/(X  -  2*RadBmax  -  2*Pspace)/ 
:       (Y  -  2*RadBmax  -  2*Pspace) 


WRITE(50,10 

WRITE(50,50 

are:  ' 

WRITE(50,40 

WRITE(50,20 

WRITE(50,20 

WRITE(50,20 

WRITE(50,20 

-  1)) 

WRITE(50,20 
WRITE(50,20 
WRITE(50,30 

-  1 
WRITE(50,30 

-  1 

WRITE(50,20 
WRITE(50,20 
WRITE(50,20 
WRITE(50,20 


'The  area  percentage  of  particles  is',  AP 
'The  sub-area  percentages  of  particles 

AF1 , AF2 , AF3 , AF4 , AF5 
'X  dimension  is',X 
'Y  dimension  is',Y 

'Min  spacing  between  particles  is',Pspace 
'Mean  particle  radius  is',  TRad/ (Pnum  -  (b 

'Total  area  occupied  by  particles',  Parea 
'Total  number  of  particles  is',  Pnum 
'Number  of  particles  along  the  border  is',b 

'Number  of  particles  in  area  interior  is', a 

'Number  particles  per  unit  area',  NPA 
'Average  nearest  neighbor  spacing  is' ,  dpavg 
'Min  nearest  neighbor  spacing  is', dpmin 
'Max  nearest  neighbor  spacing  is', dpmax 
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WRITE (50,  20)  'Max  particle  radius  is',Radmax 
WRITE(50,20) 'Min  particle  radius  is',Radmin 
PRINT*, '  The  area  percentage  of  particles  is',  AP 
PRINT' (T6,A,T45,5F4.1) ' , 'The  sub-area  percentages  of 
: particles  are ' , AF1 , AF2 , AF3 , AF4 , AF5 
PRINT' (T6,A,T45,5F6.2) ' , 'The  area  occupied  by  particles 
: in  each  sub-area  is...',Al,  A2 ,  A3,  A4,  A5 
PRINT*,'  X  dimension  is',X 
PRINT*,'  Y  dimension  is',Y 

PRINT*,'  Min  spacing  between  particles  is',Pspace 
PRINT*,'  Mean  particle  radius  is',  TRad/ (Pnum  -  (b  -  1)  ) 
PRINT*,'  Total  area  occupied  by  particles  is',  Parea 
PRINT*,'  Total  number  of  particles  is',  Pnum 
PRINT*,'  Number  of  particles  along  the  border',  b  -  1 
PRINT*,'  Number  of  particles  in  area  interior',  a  -  1 
PRINT*,'  Number  particles  per  unit  area  is',  NPA 
PRINT*,'  Average  nearest  neighbor  spacing  is',  dpavg 
PRINT*,'  Min  nearest  neighbor  spacing  is',dpmin 
PRINT*,'  Max  nearest  neighbor  spacing  is',dpmax 
PRINT*,'  Max  particle  radius  is',Radmax 
PRINT*,'  Min  particle  radius  is',Radmin 

10     FORMAT(T6,A,T50,F10.4) 

20     FORMAT(T6,A,T50,F10.4) 

30     FORMAT(T6,A,T45,I10) 

40     FORMAT (5F6. 2) 

50     FORMAT (A) 
END 
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